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I. INTRODUCTION 


Electron bombardment Ion thrusters have been proposed for both 
primary and auxiliary space propulsion applications for nearly two-and- 
a-half decades. During this time, thruster performance requirements 
have varied according to the missions of current interest. Also during 
this time, the development of thruster designs capable of meeting these 
requirements has been largely experimental. That is, thruster develop- 
ment has generally been accomplished by a procedure in which the design 
parameters that influence thruster performance are physically varied 
until an acceptable configuration is obtained. 

This procedure has several Inherent limitations. First of all, it 
can be time consuming, especially If a large number of parameters is 
involved. Secondly, an optimum configuration may not be found. That Is, 
once a configuration capable of fulfilling the mission requirements is 
identified, the procedure is generally terminated even though this may 
not result in the best configuration possible. Finally, changes in the 
missions of current interest to those characterized by different 
thruster requirements necessitate that the iterative experimental pro- 
cedure be repeated. 

Consequently, there Is a need for the development of an analytical 
model which describes the effects of thruster design variables and 
operating parameters on thruster performance. Such a model should, as 
a minimum, be capable of providing guidance for the iterative procedure 
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described above, and ideally would be capable of describing exactly how 
a thruster should be designed to achieve a given set of performance 
requirements. 

Many models of discharge chamber operation have been developed 
over the last 24 years. A brief discussion of these models is given in 
the next section. In general, the complexity of the processes taking 
place in the discharge chamber of an ion thruster together with the 
relative ease with which new thruster designs can be tested experi- 
mentally has resulted in a situation in which theoretical understanding 
has lagged experimental developments. 

The objective of this research is to improve the theoretical 
understanding of ion thruster operation by providing a simple physical 
model of the processes affecting thruster performance. Additional con- 
straints on this model are that it should be easy to use, yet general 
enough to be applicable to a wide range of thruster configurations and 
operating conditions. 


Background 

Analytical modeling of electron bombardment ion thrusters has been 
on going more or less continuously since their introduction by Kaufman 
[1,2] in 1960. Milder [3] provides a survey of modeling efforts made 
through 1969. In this survey he concluded that, "...our knowledge and 
understanding of the physics of these plasma? is far from complete." 

In addition, he concluded that the usefulness of these efforts was 
limited by either the slmplyfing assumptions required to make the 
nroblem tractable or the difficulty of applying less simplified models. 


3 


Kaufman [4] presents a discussion of thruster technology as of 
1974 In which the latest theories of thruster operation are described. 

He concluded, as had Milder earlier, that the most successful effort 
to date was the semi empirical approach proposed by Masek [5] in 1969. 

The theory of Masek Is semi empirical In that it requires as Input a 
detailed knowledge of the plasma properties inside the thruster dis- 
charge chamber. These properties are generally obtained using a 
Langmuir probe together with a data reduction procedure that is both 
tedious and until recently [6,7] of only limited accuracy. 

Since 1974 a number of analytical models have been developed, 
directed toward different aspects of ion thruster discharge chamber 
operation [7-19]. Of these, the models proposed in references 8 
through 12 are extensions of Masek's modeling technique in that they 
require detailed Langmuir probe data as inputs. References 16 through 
18 present models which do not require probe data, but which instead 
consist of a complex set of equations which must be solved iteratively 
by a computer. This complexity makes these models difficult to apply 
and limits their usefulness in providing guidelines to improved 
thruster designs. 

Discharge chamber models of ion sources for neutral beam injectors 
have also been proposed recently [20-22]. The operation of these 
sources is in many respects very similar to that of ion sources for 
space propulsion applications. This similarity was enhanced by the 
recent switch in space thruster design [23,24] to discharge chambers 
characterized by the same high magnetic flux density cusped fields used 
in neutral beam injectors. References 20 through 22 present relatively 
simple models that provide valuable insight to discharge chamber 
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operation, but also require plasma property data as inputs. Further, 
these models are not developed to the point where the performance of a 
given discharge design can be calculated directly. 

In summary, the situation at the present time Is remarkably 
similar to the way It was In 1969 as described by Milder [3]. Many 
additional models have been proposed since 1969 providing a significant 
improvement In the overall understanding of discharge chamber processes, 
however, most of these models still require detailed plasma property 
data as inputs. Those that do not are either over simplified, result- 
ing in a loss of generality and usefulness, or too complex to be 
applied easily. Thus, there is a need for the development of a model of 
ion thruster performance that is easy to use, does not require plasma 
data as inputs, and yet is general enough and accurate enough to serve 
as a guideline for the design of improved thruster configurations. 

Overview 

This dissertation is organized in the following manner. First, a 
brief review of ion thruster operation is given in Chapter II. In this 
chapter, the dominant mechanism affecting thruster performance are 
identified and discussed qualitatively. Since this report is primarily 
concerned with ion thruster discharge chamber operation, detailed dis- 
cussions of other thruster components such as the ion accelerator 
system and cathodes, etc., will not be given. These components will be 
discussed only in regard to their effect on the operation of the main 
discharge chamber. 

In Chapter III the analytical derivation of the equations com- 
prising the thruster performance model proposed in this investigation 


Is presented. This model consists of a set of algebraic equations, 
each of which describes the behavior of a different performance param- 
eter. For example, a single equation for the thruster performance 
curve (variation of beam Ion energy cost with propellant utilization 
efficiency) Is developed. Other equations are developed that allow the 
average values of the discharge chamber plasma properties to be calcu- 
lated as functions of the thruster operating point. These plasma 
properties Include: the primary electron density, the prlmary-to- 
Maxwelllan electron density ratio, the Maxwellian electron temperature, 
and the ratio of doubly-to-slngly charged Ion currents in the beam. 

In Chapter IV, the experimental apparatus and procedures used to 
investigate the validity of the proposed model are described. Experi- 
mental results, comparison of these results with the predictions of the 
model, and a discussion of this comparison are also given in this 
chapter. 

Three applications of the thruster performance model are discussed 
In Chapter V. The first of these illustrates the effect of thruster 
design and operating parameters on the standard performance curve. The 
second application describes how the model can be used to extrapolate 
data taken without ion beam extraction to obtain the performance curve 
applicable to operation with beam extraction. Finally, an illustration 
of how the model can be applied to mission analysis calculations Is 
presented. In this application, the effects of thruster design param- 
eters and propellant utilization efficiency on the deliverable payload 
fraction are discussed for an earth orbit raising mission (low earth 
orbit to geosynchronous earth orbit). 


* aaeawdut vrr.i wm 




The major conclusions of this Investigation are summarized In 
Chapter VI. SI units are used throughout this report with the excep- 
tion that energy Is frequently given In units of electron volts. In 
addition, for thermal electrons the energy quantity KT. where k is 
Boltzmann's constant and T Is the temperature will be given as eT 
where e Is the electronic charge and T„ Is the electron temperature In 



II. ION THRUSTER OPERATION 


A contemporary "ring cusp" electron bombardment ion thruster is 
shown schematically in Fig. 1. This thruster consists of a cylindrical 
steel structure bounded at one end by a circular steel plate and at the 
other end by a set of two closely spaced grids. The volume bounded by 
this structure is referred to as the discharge chamber. The two grids 
at the end of the discharge chamber have a matched matrix of holes 
through which the ion beam is extracted. The Inner grid is referred to 
as the screen grid while the outer one is called the accelerator grid. 

During thruster operation, neutral propellant gas is injected into 
the discharge chamber. The preferred propellant gas has, for the most 
part, been mercury vapor because of its large atomic mass, low ioniza- 
tion energy, and its liquid phase storability. Cesium vapor has also 
been used. However, at the present time, the propellants of most in- 
terest are the rare gases argon, krypton and xenon. In this investiga- 
tion, only these rare gases are considered. The propellant gas is 
assumed to fill the discharge chamber uniformly during thruster 
operation. 

A hollow cathode serves as the source of electrons for the dis- 
charge chamber. Hollow cathodes have replaced both refractory metal 
and oxide cathodes in ion thrusters primarily because they exhibit very 
long lifetimes and can be restarted a number of times even after 
exposure to air (an important capability if the thruster is to be 
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preflight tested). A detailed discussion of hollow cathodes Is given 
by Siegfried [25]. 

The electrons emitted by the cathode In Fig. 1 are accelerated by 
an electric field adjacent to the cathode Into the discharge chamber. 
This electric field is established by biasing the discharge chamber 
walls (with the exception of the screen grid) 30 to 50 volts positive 
of the cathode by means of an external DC power supply (called the 
discharge or anode supply). Electrons which have undergone this accel- 
eration are called primary electrons. The energy of these primary 
electrons is determined by the voltage difference applied between the 
anode and cathode. The magnitude of this voltage difference is chosen 
so that the collision cross section for ionization by the primary 
electrons is large while at the same time the collision cross section 
for the production of multiply charged ions by the primary electrons is 
small . 

A second group of electrons originates from the inelastic inter- 
action of the primaries with the neutral propellant atoms. These 
interactions reduce the energy of the primary electrons. In this in- 
vestigation, an electron is no longer considered a primary electron if 
it has had at least one inelastic collision. Ionization, the inelastic 
collision process of primary interest here, results in the release of 
low energy secondary electrons. Primary electrons, that have had their 
energy degraded by inelastic collisions, and secondary electrons 
released by ionization, thermal ize to form an electron population with a 
nearly Maxwellian energy distribution (characterized by a temperature on 
the order of a few eV). It is possible for both electron populations 
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to exist simultaneously In the discharge chamber due to the low Inter- 
action rate between the primary and Maxwellian electrons [26]. 

At typical discharge chamber neutral atom densities M0 18 cnr 3 ), 
the Ionization mean free path for primary electrons In neutral atoms Is 
on the order of meters while typical discharge chamber dimensions are 
on the order of tens of centimeters. For this reason, a magnetic field 
is employed to restrict the direct access of primary electrons to the 
anode. The magnetic field for the thruster of Fig. 1 is created by 
the rings of magnets of alternating polarity located along the back and 
sides of the discharge chamber. The configuration and strength of the 
magnetic field has a substantial effect on thruster operation and con- 
sequently has been the subject of numerous studies [6,7,13,24,27-33]. 
Contemporary thruster designs are operated with the entire discharge 
chamber housing (except for the screen grid) at anode potential. Mag- 
netic field lines at the field cusps, thus, terminate on anode potential 
surfaces. This allows electrons to be lost to the anode by travelling 
along magnetic field lines as well as by diffusing across them. To 
adequately restrict the flow of primary electrons along the field lines 
to the anode, magnetic flux densities on the order of 0.1 tesla are 
required. 

The discharge chamber magnetic field reduces the probability that 
a primary electron will be collected by the anode without first having 
had an inelastic collision with a neutral propellant atom. This prob- 
ability is a function of the thruster size, the discharge chamber mag- 
netic field configuration, the cathode location and the neutral atom 
density. As the neutral density decreases, the probability that a 


primary electron will be lost to the anode without having an inelastic 
collision increases. 

For the hypothetical case of a zero neutral atom density, this 
probability is one. That is, all primary electrons emitted by the 
cathode will be collected by the anode and none will have inelastic 
collisions. In this case, each primary electron will, on the average, 
travel a certain distance through the discharge chamber before being 
collected by the anode. This distance is a characteristic of the 
thruster geometry, magnetic field configuration and cathode location, 
and is called the primary electron containment length [34]. With this 
definition, the probability of primary electron loss to the anode may 
be expressed as a function of the ratio of the primary electron con- 
tainment length to the mean free path for primary electron-neutral atom 
inelastic collisions [34]. The loss of primary electrons to the anode 
constitutes a loss of discharge energy. Consequently, ion thruster 
performance is strongly dependent on the probability with which primary 
electrons are lost [34,35]. 

The plasma produced within the discharge chamber will typically 
assume a potential a few volts positive of the anode. Thus, a 
potential sheath will exist at all plasma boundaries. The magnitude of 
this sheath at cathode potential surfaces depends primarily on the 
magnitude of the discharge voltage applied between the cathode and 
anode. The sheath potential at cathode surfaces is sufficiently nega- 
tive to reflect all but the most energetic electrons In the tail of the 
Maxwellian distribution. Consequently, the vast majority of electrons 
In the plasma can leave the discharge chamber only at anode potential 
surfaces. 
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Of the Ions produced In the discharge chamber, some will reach the 
discharge chamber walls. Those that do, recombine with electrons there 
and return to the plasma as neutral atoms. Most of the Ions that reach 
the grid system, on the other hand, are accelerated by the electric field 
between the screen and accelerator grid to form the Ion beam. The 
fraction of the total ion current produced that is extracted Into the 
beam is called the extracted ion fraction. Because it does little good 
to produce Ions In the plasma only to have them recombine at the dis- 
charge chamber walls, It Is clearly desirable to have thruster designs 
In which the extracted Ion fraction Is as large (close to one) as pos- 
sible [34]. The electric field between the grids Is established by 
biasing the thruster body positive of ground potential (on the order of 
1000 volts) and biasing the accelerator grid several hundred volts nega- 
tive of ground potential. The final velocity of the Ions In the beam 
is determined by the sum of the positive potential applied to the 
thruster cathode and the discharge voltage. 

Electrons are Injected into the positive Ion beam by a cathode 
(called a neutralizer) positioned downstream of the accelerator grid In 
order to space charge and current neutralize the beam. The negative 
accelerator grid prevents these electrons from backstreamlng to the 
positive thruster body. 

During thruster operation, the propellant gas In the discharge 
chamber is only partially Ionized, with the Ion density typ ically an 
order of magnitude smaller than the neutral atom density. Thus, one 
might expect the neutral flux through the grids to be greater than the 
ion flux. This, however. Is not the case since the neutral flux Is 
determined by free-molecular flow toward the grids at a temperature 
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governed by the mean discharge chamber wall temperature [36] 400K). 

The Ions, on the other hand, tend to move toward the grid system at the 
Bohm [37] velocity, which Is governed by the much higher electron 
temperature (typically 5 eV or ^ 58,000K). Thus, the Ion flux through 
the grids is usually greater than the neutral flux even though the loti 
density is substantially smaller than the neutral density. 

The ratio of the beam current (assuming singly charged ions only) 
to the total propellant flow rate (expressed in units of equivalent 
amperes) is known as the propellant utilization efficiency. For a 
constant total propellant flow rate, increasing the propellant utiliza- 
tion efficiency generally requires an increase in the average energy 
required to produce a beam ion. This variation of the beam ion energy 
cost with the propellant utilization is known as a performance curve. 

It will be shown in the chapters that follow that the dominant mechan- 
isms affecting the performance curve are the loss of primary electrons 
to the anode and the extracted ion fraction. 


III. THEORETICAL DEVELOPMENT 


Thruster Performance Model 

In this section, a model is developed which results directly in a 
single algebraic equation for the variation of the beam Ion energy cost 
with the propellant utilization efficiency for high flux density cusped 
magnetic field thrusters. This is made possible by formulating the 
model in terms of the average energy required to produce an Ion in the 
discharge chamber plasma and the fraction of ions produced which are 
extracted into the beam. A key feature of the model is that it 
provides a simple technique for the calculation of the average dis- 
charge plasma ion energy cost as a function of the propellant flow rate 
and propellant utilization efficiency. 

Assumptions and Limitations 

The development of the model assumes steady state discharge 
chamber operation. The model is applicable to discharge chambers which 
produce low pressure, partially ionized, optically thin plasmas In 
which neutral densities are In the range of 10 18 to 10 19 m' 3 and plasma 
densities range from 10^ to lO* 7 nf 3 . In addition, Maxwellian elec- 
tron temperatures and primary electron energies should range from 
1 to 10 eV and 30 to 50 eV, respectively. 

Electron energy losses due to elastic collisions with ions or 
neutral atoms are neglected. This can be justified because the average 
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energy loss per encounter Is proportional to the ratio of masses* and 
this ratio is small lO" 5 ). Electron energy losses resulting from 
Inelastic collisions with ions are also neglected because of the low 
ion density relative to the neutral density. This should not cause 
a significant error, except perhaps at very high propellant utiliza- 
tions. Only singly charged ions are assumed to be produced, and the 
effect of metastable atomic states on ion production is neglected. 

Primary electron thermal izati on resulting from collisions with the 
background Maxwellian electrons also is neglected. Primary electron 
behavior is assumed to be limited to either Inelastic collisions with 
neutral atoms or direct loss to anode potential surfaces. A primary 
electron Is considered to join the Maxwellian electron population after 
having one Inelastic collision. 

Electrons are assumed to be constrained by the plasma sheaths so 
they are able to leave the plasma only at anode potential surfaces. 

Ions and photons (emitted by the de-excitation of excited propellant 
atoms) are assumed to be lost across all plasma boundaries. The 
assumption of a low pressure discharge implies that ion-electron re- 
combination should be wall controlled, consequently, volume ion recombi- 
nation is neglected. 

The neutral atom density is assumed to be uniform throughout the 
discharge chamber, and free molecular flow is assumed to apply to the 
neutral atom flux through the accelerator grid system. 

Beam Ion Energy Cost 

The average energy cost per beam ion Is defined as, 

e B 5 ( J D " V V J B * C eV/beam 1on 3 (1) 


* 
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where all symbols are defined in Appendix C. The beam current (Jg) In 
Eq. 1 Is subtracted from the discharge current (Jp) so that the energy 
that goes into accelerating the beam ions through the discharge voltage 
is not charged to the beam ion energy cost. The numerator on the 
right-hand-side of Eq. 1 represents the power used to operate the dis- 
charge chamber. 

In a similar manner, the average energy expended in creating ions 
in the discharge chamber plasma may be defined as, 


e p 3 £0 d “ ( j b + J C>3 V J P * CeV/plasma ion]. (2) 


By analogy to Eq. 1, the "Jg + Jp" term is subtracted from the dis- 
charge current so that the energy that goes into accelerating these 
ions out of the discharge chamber plasma into the chamber walls or the 
beam is not included in the plasma ion production cost. Rearranging 
Eq. 2 yields, 


= ( J D~ J b) V D / \ J C V D 
P \ 


(3) 


For steady state operation, the total ion current produced (J p ) must be 
equal to the total ion current leaving the plasma. For ion thruster 
discharge chambers, ions can only leave the plasma by going to cathode 
potential surfaces, anode potential surfaces or by being extracted into 
the beam. Thus, the total ion current produced is given by*, 


J P a J B + J C + J A 

Dividing this equation through by J p yields, 


( 4 ) 


* Equation 4 is a statement of continuity for the ions. 


T 
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1 ■ U + f C + f A ’ 


(5) 


where. 


f B 5 J B /Jp , f c 5 J C /J P and f A 5 J A /J P 


( 6 ) 


The fractions f B . f c and are, in order, the extracted ion fraction, 
the fraction of ion current produced that goes to cathode potential 
surfaces, and the fraction of ion current produced that goes to anode 
potential surfaces. Using the definitions of f„ and f c in Eq. 3 along 

with Eq. 1 yields, 

f f V (7) 

£ p * e B f B ‘ f C V D * 

Solving this equation for e B gives the result, 


e B “ e p /f B + f C V D /f B ‘ 


( 8 ) 


This equation describes the be.. Ion energy cost as a function of the 
plasm ion energy cost U p >, the extracted Ion fraction (f B >. «e 
fraction of Ion current to cathode potential surfaces (f c ) and the 

charge voltage (Vq). 

The first tern on -he right-hand-side of Eq. 8 represents the 
energy loss associated with producing Ions In the discharge chanter and 
extracting only a fraction of the. Into the beam. Ions which are not 
extracted Into the bea. go to the walls of the discharge chanter where 
they recombine. The resulting atoms must then be re-ionized before 
they can contribute to the bea. current. The factor l/f„ may »e 
interpreted as the average number of times that a beam Ion undergoes 
ionization from a neutral state before being extracted Into the beam. 


rz 
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The second term on the right-hand-side of Eq. 8 represents the 
energy wasted in accelerating plasma ions into Interior cathode po- 
tential surfaces. This process is undesirable because it results in 
both a discharge energy loss and in the sputter erosion of these 
surfaces. 

To generate performance curves using Eq. 8, one must be able to 
specify the behavior of each of the terms on the right-hand-side of 
this equation as a function of the propellant utilization efficiency. 

Plasma Ion Energy Cost 

The plasma ion energy cost parameter (e p ) appearing in Eq. 8 and 
defined by Eq. 2 reflects all mechanisms of energy loss from the dis- 
charge chamber except for the acceleration of ions out of the plasma 
through the discharge voltage. Specifically, e p Includes energy losses 
due to the following phenomena: direct primary electron loss to the 
anode, Maxwellian electron collection by the anode, excitation of 
neutral atoms, excitation of ionic states (which will be neglected) and 
hollow cathode operation. To derive an expression for the plasma ion 
energy cost as a function of the propellant utilization, a power 
balance is made on the discharge chamber plasma represented in Fig. 2. 
The boundary of the control volume for this power balance is defined by 
the plasma sheath edge. The primary electrons are assumed to be accel- 
erated from a cathode region plasma potential that is V^. volts above 
cathode potential to the potential of the bulk plasma which is assumed 
to he near that of the anode. In addition, if it Is assumed that only 
the discharge power supply is used to sustain the discharge, then the 
rate at which energy is supplied to the discharge chamber plasma by the 
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Figure 2. Disharge Plasma Power Balance Schematic 
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primary electrons Is given by J E (V p - V c ). The "missing" power J E V C Is 
used to operate the hollow cathode. 

Energy Is lost from the plasma primarily by the flux of four types 
of energy carriers across the plasma boundaries: Ions, photons (emitted 
by de-excitation of excited propel lam: atoms), Maxwellian electrons, 
and primary electrons. The Ions and photons are lost to all Interior 
thruster surfaces whereas the Maxwellian and primary electrons are 
assumed to be lost to the anode surfaces only. In steady state, the 
rate of energy supplied to the plasma must be equal to the rate at 
which it Is lost, thus. 


J £ (V 0 -V c ) - JpU t + j JjUj + J M e M ♦ J L (V 0 -V C ) , (9) 

where the suimatlon Is over the set of excited neutral states. 

Dividing Eq. 9 by the ion production current (u p ) and recognizing that 
the emission current (J^) related to the discharge current by *, 


J E = J D ‘ ( J C + * 


(10) 


allows Eq. 9 to be written as. 
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. °M e M + J L V 0 J F V C J L V C 
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where Eq. 2 has been used. The rate at which the jth excited state Is 
produced (expressed as a current) Is given by, 


J j - en o"e <s jV * • 


( 12 ) 


* Equation 10 Is a statement of continuity for the electrons. 





where «jjV e > represents the product of jth excitation collision cross 
section and the electron velocity averaged over the entire electron 
speed distribution - including the primary electrons. Similarly, the 
ion production current Is given by. 
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Substituting Eqs. 12 and 13 Into Eq. 11 yields, 
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where the parameter e Q is defined by. 
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This parameter accounts for the energy that is expended in ionization 
and excitation reactions and may be calculated in the manner described 
by Dugan and Sovle [38]. The third term on the right-hand-side of 
Eq. 14 may be written as. 


J L V D _ J L 


J E V D 


(16) 


u p u p P 

where the last step was made using Eqs. 2 and 10. The ratio J L /J is 
simply the fraction of primary electrons emitted by the cathode which 
are collected by the anode before having any inelastic collisions. This 
fraction may be given by the survival equation [39] as, 


J. -ant 
L . _ °o o e 

¥ ' 6 


(17) 


where is the average distance a primary electron would travel in 
the discharge chamber before being collected by the anode - assuming it 
had no inelastic collisions (i.e. the primary electron containment 
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length), Is the total Inelastic collision cross section for 
electrons at the primary electron energy and n Q Is the neutral atom 
density. Combining Eqs. 14, 16 and 17 yields. 







o e 


(18) 


The current of Maxwellian electrons to the anode may be given as the 
sum of the secondary electrons liberated in the ionization process and 
the thermal ized primary electrons, thus. 


= + = + * (19) 


Using Eqs. 17 and 19 in 18 yields, 


-a' n £ e u -a'n i, 

E P ' e 0 + [ J p ♦ >lr + ' P t« 006 


-a' V- 
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( 20 ) 


Finally, using Eqs. 2 and 10 and solving Eq. 20 for £p results in. 



( 21 ) 


o o e 

The factor 1-e may be interpreted as the probability that a 

primary electron will have an inelastic collision before being collected 
by the anode. This is analogous to the "fast neutron non-leakage prob- 
ability" used in nuclear reactor physics [40]. 

The neutral density (n Q ) may be expressed in terms of the pro- 
pellant flow rate and propellant utilization by equating the rate at 
which propellant enters and leaves the discharge chamber, i.e., 


m = J B + 




( 22 ) 


where m and n are In units of equivalent amperes*. The neutral flow 
rate from the thruster may be expressed using the theory for free 
molecular flow through a sharp-edged orfice [39], 

% = * n o ev o A g *o ‘ (23) 

Combining Eqs. 22 and 23 yields. 
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where the propellant utilization defined by n u 
Eq. 21 may be written as. 


- Jg/ii) was used. 


(24) 

Thus, 


where. 
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and. 
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£p * ■ < 27 > 
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Equation 25 Is a very simple relationship which can be used to calculate 
the plasma Ion energy cost as a function of the propellant utilization. 
Experimental results, which will be presented In Chapter IV, Indicate 
that under many conditions the parameters C„ and s‘ may be taken to 
be Independent of the propellant utilization. 

Substitution of Eq. 25 Into Eq. 8 yields the following single 
equation describin g the performance of a given thruster design, 

+ Equation 22 is a statement of continuity for the propellant. 
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For design purposes, the parameters fg and fg in addition to C and 
e p “ay be taken to be independent of the utilization and flow rate. 
These parameters do, however, depend strongly on the thruster design. 
Indeed, these four parameters determine the performance of a given 
thruster design. 

The quantity C Q reflects the degree to which primary electrons 
interact with neutral atoms. Thus, it is referred to as the primary 
electron utilization factor. This factor depends on the quality of the 
primary electron containment (through * e ), the quality of the contain- 
ment of neutral atoms (through A g , <t» 0 and v Q ) and the propellant gas 
properties (through a Q and v Q ). Recall that the primary electron con- 
tainment length may be interpreted as the average distance a primary 
electron would travel in the discharge chamber before being lost to the 
anode— assuming it had no inelastic collisions. Magnetic fields in all 
discharge chamber designs serve the function of Increasing this length. 
Although an effective means of determining * remains to be developed, 
it is believed that this parameter is a function primarily of the 
thruster geometry, magnetic field configuration and cathode location. 

Equation 25 suggests that the plasma ion energy cost should, 
through the factor C Q , be a function of: 

1. The propellant, which determines the total Inelastic collision 
cross section (o^) and atomic mass (which affects the thermal neutral 
velocity, v Q ). 

2. The wall temperature, which affects the neutral velocity. 
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3. The transparency of the grids to neutrals U 0 ). 

4. The area through which the beam Is extracted ( A g ). 

5 . The discharge voltage, which determines the primary electron 
energy and. thus, effects the value of the cross section (o 0 ). 

The baseline plasma Ion energy cost U p ). defined by Eq. 27, 
depends on a number of energy loss mechanisms Including: the relative 

amount of energy expended In excitations compared to Ionization of 
neutral atoms through V the average energy of the Maxwellian electrons 
which leave the plu. U„) and the efficiency with which the hollow 

cathode operates. 

The cathode efficiency Is reflected In the value of \ which 
represents the plasma potential from which the electrons are supplied. 
Inefficient cathode operation results In high values of V c and corres- 
pondingly poor overall thruster performance. For thermionic cathodes 
v c = 0, however, additional heater power must be supplied to effect Us 
operation. For thrusters equipped with a cathode pole piece/baffle 
assembly. V c should be taken as the plasma potential In the cathode 
discharge region (l.e., the region between the cathode and the baffle). 

In this case, the power represented by J E V C goes Into both the operation 
of the hollow cathode and the operation of the cathode region discharge. 
The resulting high values of V c , In this case, would be expected to 
produce poorer overall thruster performance. Elimination of the 
separate cathode discharge region should, therefore, Improve the per- 

formance. 

The parameter c o , defined by Eq. 15, in general Is a function of 
the entire electron energy distribution Including the primary electrons. 
The parameter cj (which contains V see Eq. 27) Is, thus, also a 
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function of the electron energy distribution. For plasmas characterized 
by a Maxwellian plus monoenergetlc electron energy distribution, a 
simple method for the calculation of eJ may be derived. Sample calcu- 
lations using this method, along with the appropriate collision cross 
section data, are given in Appendix A. The experimental results given 
in Chapter IV indicate that, under many conditions, the baseline plasma 
ion energy cost (e*) may be taken to be a constant. The calculations 
given in Appendix A demonstrate that this experimental observation is 
predicted by the model. 

Although the model cannot yet be used to predict the performance 
of completely new thruster designs, it provides a clear physical pic- 
ture of the phenomena affecting the performance. Equation 25 describes 
the plasrr ion energy cost in terms of the loss of primary electrons to 
the anode. At high values of the neutral density parameter, the neutral 
density in the discharge chamber is large, and the probability is high 
that all the primary electrons will undergo Inelastic collisions with 
neutral atoms and none will be lost directly to the anode. In this 
case, the discharge chamber will be producing ions for the minimum or 
baseline plasma ion energy cost, e*. 

As the beam current is increased (at a constant propellant flow 
rate), the propellant utilization increases causing the neutral density 
parameter (and thus the neutral density) to decrease (see Eq. 24). The 
decrease in neutral density Increases the likelihood of a primary elec- 
tron reaching the anode without having an Inelastic collision. This 
direct loss of primary electron energy increases the overall plasma ion 
energy cost according to Eq. 25 and consequently increases the beam ion 
energy cost according to Eq. 28. The shape of the performance curve is 
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largely determined by this direct loss of primary electrons. Thus, any 
design change which decreases the likelihood of the direct loss of 
primary electrons (without decreasing the extracted ion fraction) 
should Improve the thruster's performance. The probability of direct 
primary electron loss is determined through the parameter C Q . This 
parameter may be increased by either increasing *, , which makes it 

c 

harder for primary electrons to escape the plasma, or by making it 
harder for neutral atoms to escape the discharge chamber. 

Calculation of Plasma Properties 

The following analysis provides a set of very simple algebraic 
equations which allow one to calculate the values of the following dis- 
charge chamber plasma properties: the average primary electron density, 

the average primary- to-total electron density ratio, the average Max- 
wellian electron temperature and the average ratio of the doubly-to- 
singly charged ion current in the beam. Each of these quantities may be 
calculated as a function of the propellant flow rate and propellant 
utilization using only information that would normally be available in 
the thruster design phase. 

Primary Electron Density 

This analysis is based on the recognition that all of the energy 
supplied to the discharge chamber plasma is supplied by the primary 
electrons. Thus, correctly accounting for the behavior of the primary 
electrons is essential in determining the average plasma properties. 

It is assumed that a primary electron can do only one of two things. 

It can either have an inelastic collision with a neutral propellant 
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atom or It can be lost directly to the anode. For the case where the 
primary electron has an Inelastic collision with a neutral atom. It 
will produce either an Ion or an excited neutral state. After such an 
inelastic collision, the energy of the primary electron is degraded and 
it is assumed that it is subsequently thermalized into the Maxwellian 
electron population. The rate at which primary electrons have in- 
elastic collisions with neutral atoms is given by the difference be- 
tween the rate at which they are supplied by the cathode (d E ) and the 
rate at which they are lost directly to the anode (J L ), i.e.. 


o E - j l 


J p ♦ 


(29) 


where Jp is the ion current produced by primary electrons and J* is 

©X 

the primary electron induced production rate of excited neutral states 
expressed as a current. The ion current produced by primary electrons 
is given by, 

J P = n o n p ea + * v p » (30) 


where o.J. is the ionization collision cross section at the primary elec- 
tron energy. Similarly, the rate of production of excited state atoms 
induced by primaries is given by, 


J 



n n ev_ V- 
o p p 


(31) 


where oj is the collision cross section for the jth excited state at 
the primary electron energy. The fraction of the primary electron- 
neutral atom inelastic collisions that produce ions is given by, 
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or, 
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J p + J e* * | °j 


( 33 ) 


Note that the sum o! + J o! Is just the total inelastic collision cross 
section for primary electron-neutral atom collisions a Q . Therefore, 
Eq. 33 may be written as. 
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(34) 


Multiplying Eq. 34 by the rate at which primary electron-neutral atom 
collisions occur (Eq. 29) yields the rate at which ions are produced by 
primary electrons, i.e., 

I 

O , 


Jp = (Jf “ j l) , 


(35) 


The term in parentheses in Eq. 35 may be written as. 


J E - J^ 88 Jj: (1 - J^/Jp) • 


(36) 


The quantity J,_/J E is the fraction of the input primary electron current 
lost directly to the anode and is given by the survival equation (Eq. 17 
written in a slightly different form), 






(37) 


Combining Eqs. 2 and 10 yields the following expression for the cathode 
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Combining Eqs. 25, 35, 36, 37 and 38 yields the following expression 
for the ion current produced by primary electrons. 
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Dividing both sides of this equation by Jp yields an expression for 
the ratio of ion current produced by primaries to the total ion current 
produced, 
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Jn 


£ P CT + 

Vo 


(40) 


Since the ratio Jp/Jp cannot be greater than one, Eq. 40 provides a 
theoretical limit for the maximum value of £p , i.e., 
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The total ion current produced may be expressed in terms of the beam 
current by using the definition of the extracted ion fraction (Eq. 6), 


3p “ Jg/^g • (42) 

In addition, using the definition of the propellant utilization the 
beam current may be written as, 

0 B 8 mn u . (43) 


Combining Eqs. 39, 42 and 43 yields, 
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Finally, equating Eq. 30 to Eq. 44 and solving for the primary electron 
density (n p ) yields, 
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(45) 


where Eq. 24 was used for the neutral density. Equation 45 provides an 
expression for the average primary electron density as a function of 
the propellant utilization. Remarkably, this expression does not 
depend on either the propellant mass flow rate or the primary electron 
utilization factor, both of which cancelled out in the analysis. The 
term in the square brackets should be roughly a constant for a given 
thruster design, propellant gas and discharge voltage. A reasonable 
estimate for each of the terms in the square bracket should be possible 
for a thruster being designed (assuming a method is developed for the 
determination of the extracted ion fraction). 


Primary-to-Total Electron Density Ratio 

Assuming quasi -neutral ity, the average ion density (n^ is equal 

to the average total electron density (n, = n + nj. In addition, the 

1 p M 

average ion density is related to the beam current by [41]. 


J 8 * 0.6n ( ev b A g * j 


(46) 


Solving this equation for n^ yields, 
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(47) 


where ^ is the transparency of the screen grid to ions and v fa is the 
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Bohm velocity. Dividing Eq. 45 by 47 and using Eq. 43 yields the 
following expression for the average ratio of primary- to- total electron 


density* 
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This equation indicates that the average primary-to-total electron 
density ratio should be a function of the neutral density parameter* 
m(l-n u ). As was the case for Eq. 45, the combination of parameters In 
the square brackets of Eq. 48 should be roughly a constant for a given 
thruster design, propellant and discharge voltage. Equation 48 also 
indicates that, all else being equal, increasing the extracted Ion 
fraction should decrease the primary-to-total electron density ratio. 


Primary-to-Maxwellian Electro n Density Ratio 

Once the primary-to-total electron density ratio is calculated, 
using Eq. 48 the primary-to-Maxwellian electron density ratio (np/n^) 
may be calculated from the following equation, 

^ , Vi > 1— . (49) 

n M 1_n p /n 1 n 1 /n p " 1 


Maxwellian Electron Temperature 

The total ion current produced (Op) is the sum of the ion current 
produced by primary electrons (Jp) and that produced by Maxwellian 
electrons (J p>M ) . 
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The Ion current produced by primary electrons Is given by Eq. 30 and 
that produced by Maxwellian electrons is given by. 


J P,M * W <0 +Vm * > 


(51) 


where the quantity < ^ + v e >m represents the product of the electron 
velocity and ionization collision cross section averaged over the 
Maxwellian electron energy distribution. This quantity is the Max- 
wellian electron rate factor for the production of ion 1 - and is given 


the symbol Q q , i.e., 

q o' < 0 .Vm • 

Combining Eqs. 50 to 52 yields. 
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Using Eqs. 42 and 43 In Eq. 53 and solving for Q* gives, 
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Substituting for the neutral density using Eq. 24 yields, 
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The Maxwellian electron density, appearing In Eq. 55, may be written as, 

n M * ( % ) n P • (56) 


Substituting this into Eq. 55 yields, 
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Equation 45 provides an expression for the primary electron density. 
Using this In Eq. 57 yields, 
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The prlmary-to-Maxwel 1 Ian electron density ratio may be found using 
Eqs. 48 and 49. Substituting these equations Into Eq. 58 yields the 

following expression for the average Maxwellian electron rate factor 
for ion production, 
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Equation 59 provides an expression for the Maxwellian electron ioniza- 
tion rate factor as a function of the neutral density parameter, 
m(l-n u ) . Once this rate factor has been calculated, It can be compared 
to a tabulation of rate factors versus electron temperature for the 
given propellant to determine the appropriate Maxwellian electron 
temperature. Equation 59 Indicates that the average Maxwellian elec- 
tron temperature should Increase with decreasing values of the neutral 
density parameter. 


Double Ion Production 

In the previous analysis, the production of doubly charged Ions was 
neglected. In this section, a simple formulation for the ratio of 
doubly-to-singly charged Ion current In the beam Is developed. The pro> 
ductlon rate of singly charged Ions expressed as a current Is given by, 


, 4+, 
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J p ■ (n oVo + " 0 % 0> e * • (60) 

where Q* Is the Maxwellian rate factor for the production of single 
ions from ground state neutral atoms and P* Is the primary electron 
rate factor also for the production of single Ions from ground state 
neutrals. Similarly, the production rate of doubly charged Ions 
expressed as a current is given by, 

v ■ ( n 0 n M r * v P c * vm r + vpC) 2e * ■ (si ) 

where n + is the singly charged ion density, and P* + are the Max- 
wellian and primary electron rate factors for double ion production 
from ground state neutrals, respectively, and and P* + are the cor- 
responding rate factors for double ion production from ground state ions. 
Dividing Eq. 61 by 60 yields. 



Assuming that the doubly charged extracted ion fraction (fj + ) is equal 
to the singly charged extracted ion fraction (fj) implies that, 
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where Jg / Jg is the average ratio of doubly-to-singly charged ion cur- 
rent in the beam. Combining Eqs. 62 and 63 yields, 
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The singly charged ion density (n + ) may be conservatively approximated 
by assuming the beam current is made up entirely of singly charged ions, 
thus, 

"♦ ' "1 = oTfi^, ' (65) 

In addition, using Eq. 43 for the beam current and Eq. 24 for the 
neutral density allows Eq. 64 to be written as. 



The first term on the right-hand-side of Eq. 66 represents the produc- 
tion of double ions from ground state neutral atoms. This term depends 
on the thruster geometry only through the primary-to-Maxwellian elec- 
tron density ratio (*ip/n M ) which is given by Eqs. 48 and 49 and the 
Maxwellian electron temperature (Eq. 59). The second term in Eq. 66 
represents double ion production from singly charged ions. This term 
is strongly dependent on the propellant utilization. The rate factors 
Q* + and Q* + may be determined once the rate factor Q* is calculated 
from Eq. 59. 
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IV. EXPERIMENTAL PROCEDURES AND RESULTS 


Apparatus 

For this investigation, the ion sources shown schematically in 
Figs. 3a and b were designed and built. Each of these sources normally 
produces a 12 cm dia. ion beam and provides the capability for measuring 
the distribution of ion currents to the beam, screen grid and internal 
thruster surfaces (with the exception of the anode). 

The magnetic field for the experimental ion source in Fig. 3a is 

established through the use of an electromagnet located on the upstream 

centerline of the discharge chamber and a number of 1.9 cm x 1.3 cm x 

0.5 cm samarium cobalt permanent magnets. These permanent magnets are 
arranged end-to-end to form ring magnets of alternate polarity in the 
manner suggested by Fig. 3a. The flux density at the surface of the 
magnets is 0.27T and the magnets are attached to the steel discharge 
chamber housing by their own magnetic attraction. This arrangement 
allows the ion source magnetic field configuration to be altered quickly 
and easily by simply adding, removing or changing the position of the 
magnets. Although many different configurations were tested, the results 
obtained were all similar, thus, only those obtained using the configura- 
tions shown In Figs. 3a and 3b will be presented. For these configura- 
tions, the upstream magnet ring is covered with a strip of 0.13 mm thick 
steel Insulated from the magnets themselves by a strip of 0.25 mm thick 
flexible mica. This is done so that the surface of this strip can be 
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maintained at anode potential while the rest of the thruster body is 
biased negative of cathode potential. The downstream magnet ring is 
uncovered. The magnetic flux density at the surface of the electro- 
magnet for the configuration of Fig. 3a can be adjusted from zero to 
approximately 0.2 T by adjusting the magnet current from zero to 
124 A. 

The main discharge chamber cathodes, for both configurations, con - 
si st of seven 0.25cmdia. tungsten wires connected in parallel and support- 
ed by two support posts that are electrically isolated from the thruster 
body. Each cathode wire is approximately 2.8 cm long so the total 
cathode length exposed to the plasma is about 19.6 cm. These seven 
short wires in parallel are used to minimize the voltage drop across 
the cathode. A voltage drop less than 3 v at the maximum heater cur- 
rent was achieved with this system. The small voltage drop across the 
cathode results in a primary electron energy distribution that more 
closely resembes the monoenergetic distribution produced by a hollow 
cathode. The cathode wires were heated using direct currents in the 
range 6 to 8 A per wire. Tests on the configuration of Fig. 3a, were 
conducted using argon, krypton and xenon propellants. Discharge volt- 
ages were varied from 30 to 50 v for argon and 20 to 40 v for krypton 
and xenon. The discharge current was adjusted through the range of 
0.5 to 5 A by controlling the heating current through the refractory 
cathode wires. 

Two ion accelerator systems were used in this study. Tne first 
accelerator system consisted of a set of dished small hole accelerator 
grids (SHAG) v»l th a cold grid separation of 0.75 mm and screen and 
accelerator grid physical open area fractions of 0.68 and 0.30, 
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respectively. The second system consisted of a set of dished large 
hole accelerator grids (LHA6) with a cold grid separation of 0.75 mm, 
and screen and accelerator grid physical open area fractions of 0.68 
and 0.57, respectively. Both accelerator systems were normally masked 
to produce a 12 cm diameter beam. One series of test was conducted, 
however, with the SHAG set masked to produce a 6 cm diameter beam. For 
the 12 cm dia. beam tests, flow rates for both argon and krypton were 
varied from 500 to 1500 mA eq. and for xenon from 250 to 1000 mA eq. 

For the 6 cm dia. beam test, the flow rates were varied from 125 to 500 
mA eq. 

The configuration of Fig. 3b is similar to that of Fig. 3a except 
that the electromagnet assembly has been replaced by a permanent magnet 
attached to the steel backplate. In addition, the source of Fig. 3b is 
2 cm shorter than the one in Fig. 3a and the upstream magnet ring is 
positioned only 2.5 cm from the downstream end rather than 6 cm. 

This source was equipped with two Langmuir probe assemblies. The 
first probe consisted of a 0.76 mm dia. tantalum wire, 4.32 mm long, 
supported from a quartz tube insulator. This probe was positioned 
along the thruster centerline approximately half way between the cathode 
assembly and the screen grid as suggested in Fig. 3b. The second probe 
was a square piece of steel, 1 cm on a side and 0.127 mm thick that was 
positioned on the surface of the upstream magnetic ring. This probe 
was insulated from the magnet ring with a piece of 0.13 mm thick 
flexible mica. 

Measurements of the doubly and singly charged ion beam components 
were made on the thruster centerline using a collimating I x B momentum 
analyzer. Details of the use of this probe are given elsewhere [42]. 
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Tests were conducted on the configuration of Fig. 3b using both 
argon and xenon propellants. Argon propellant flow rates were varied 
over the range 350 to 1500 mA eq. at discharge voltages of 30 to 50 v. 
Xenon propellant flow rates were varied over the range 250 to 1000 mA 
eq. at discharge voltages of 20 to 40 v. Discharge currents for both 
propellants were varied over 0.5 to 4.0 A. 

All tests were conducted in a 1.2 m dia. x 4.6 m long vacuum test 
facility. Indicated tank pressures ranged from % 2 x 10“ 6 Torr with 
no flow to 3 x 10“ 5 Torr at a flow rate of 1500 mA eq. of krypton. 

Procedure 

The following set of experiments, conducted using the thruster 
configuration of Fig. 3a, was designed to test the suitability of the 
model developed in Chapter III to predict the functional dependence of 
the plasma ion energy cost, e p , on the neutral density parameter, 
m(l-n u ) • The model predicts that the plasma ion energy cost should 
behave according to Eq. 25 with the primary electron utilization factor 
(C Q ) and the baseline plasma ion energy cost (e*) given by Eqs. 26 and 
27, respectively. The value of E p may be determined experimentally 
through the use of Eq. 2, which is repeated here for convenience, 

+ J B^ V D 

£p ' ^ 

Note that the power used to operate the thermionic cathodes is not in- 
cluded in Eq. 2. In order to use Eq. 2, one must be able to measure 
each of the parameters on the right-hand-side of the equation. Measure- 
ment of the discharge current, discharge voltage and beam current is 
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straight forward, and was accomplished using the experimental Instru- 
mentation described In Appendix B. 

To measure the Ion currents Jq and J p , the thruster configuration 
of Fig. 3a was operated with only the upstream magnet ring at anode 
potential. All other interior discharge chamber surfaces (with the 
exception of the cathode support posts) were biased approximately 30 v 
negative of cathode potential to repel the discharge chamber electrons. 
At this bias, the current to these surfaces consists only of the incom- 
ing ion current. If the ion current to the cathode support posts and 
cathode wires is neglected, then the ion current measured in the manner 
described above is equal to J^. 

To determine J p , it is noted that the total ion production current 
is the sum of the ion currents leaving the plasma as given by Eq. 4. 
Since the anode is exposed to the plasma only at a magnetic field cusp, 
the effective area for ion loss to this surface is expected to be less 
than the physical area [43-46]. Rough calculations indicate that the 
ion current to the anode with this configuration should be less than a 
few precent of the total production current. Thus, J p may be approxi- 
mated as the sum of the ion currents to the beam (including the impinge- 
ment current) and to the negatively biased discharge chamber surfaces 
(including the screen grid). 

Complete sets of data consisting of the beam current, propellant 
flow rate, propellant utilization and total ion production current were 
collected over the range of operating conditions discussed earlier with 
the electromagnetic current held constant at 57 A for the thruster of 
Fig. 3a. At each condition tested, the thruster was operated at flow 
rates of 500, 750, 1000 and 1500 mA eq. for argon and krypton 
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propellants, and 250, 500, 750 and 1000 mA eq. for xenon. At each flow 
rate, the discharge voltage was held constant while the discharge cur- 
rent was varied. Increasing the discharge current, increases the beam 
current and propellant utilization, thus causing the neutral density 
parameter m(l-n u ) to decrease. By operating in this manner, the full 
range of neutral density parameters from close to zero to nearly 1500 
mA eq. could be investigated. Thus, the plasma ion energy cost was de- 
termined from Eq. 2 for a wide range of operating conditions. The ex- 
tracted ion fraction, fg, defined by Eq. 6, was also computed from 
measured ion currents over this same range of conditions. 

A second set of experiments was performed using the thruster 
configuration of Fig. 3b. These tests were designed to investigate the 
ability of the model to predict the behavior of various plasma prop- 
erties with variations in thruster operating conditions. The procedure 
for these tests was the same as for the previous set with the addition 
that Langmuir probe measurements of the plasma properties on the thruster 
centerline and at the anode surface were made at each operating point 
tested. The ratio of doubly to singly charged ion beam currents along 
the thruster centerline were also measured at each operating condition, 
using E x B probe. 


Experimental Results 

Plasma Ion Energy Cost 

Measurement of the plasma ion energy cost, for operation of the 
thruster configuration of Fig. 3a with argon propellant at a 50 v 
discharge over the range of neutral flow rates from 500 to 1500 mA eq., 
yielded the results shown in Fig. 4a. Here, the measured values of the 
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plasma ion energy cost (e p ) are plotted as a function of the neutral 
« 

density parameter m(l-n u ), as suggested by Eq. 25. The solid line in 
Fig. 4a is the curve given by Eq. 25 when the parameters C Q and e* have 
been selected to give the best fit to the data. The parameter e* was 
taken to be the value of e p measured at large values of the neutral 
density parameter. The justification for this selection can be under- 
stood by considering Eq. 25, which shows that when C Q m(l-n u ) is large, 
the exponential term is small compared to unity and one obtains 
e p = e p . Having established the value of the baseline plasma ion 
energy cost, the value of the primary electron utilization factor (C Q ) 
is varied until the best fit Is obtained. The agreement between the 
functional form of Eq. 25 and the experimental data is seen to be quite 
good. This indicates that the parameters C Q and e* may be taken to be 
Independent of the neutral density parameter. 

A value of C Q = 3.1 (A eq.)" 1 which is applicable to the ion source 
of Fig. 3a operating at the conditions defined in the legend for Fig. 4a 
has now been established. New values of primary electron utilization 
factor, applicable to other operating conditions, may now be calculated 
from this value using Eq. 26. For example, changing grid sets from SHAG 
to LHAG should change the value of C Q through the parameter <f> Q which is 
the effective transparency of the grids to neutral atoms. This effec- 
tive transparency parameter may be calculated for each grid set accord- 
ing to the equation, 


*o 


s* 
s + * 


( 66 ) 


where 4> s and <j> a are the modified transparencies for the screen and 
accelerator grids, respectively. These modified transparencies may be 
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calculated as the physical open area fraction of a grid times the 
appropriate Clausing factor [47], For the two grid sets used in this 
study, 

^o^SHAG * °* 16, ^o^LHAG = 0,27 * ( 67 ) 

Thus, the new value of the primary electron utilization factor appli- 
cable to the LHAG optics with all other conditions held constant is 
given by. 


^o^new 


^o^SHAG 

( Vlhag 


^o^Fig. 4a 


( 68 ) 


which yields (C Q ) new= 1>8 (A eq .)-i 

The measured values of e p , obtained under the same set of condi- 
tions defined in the legend of Fig. 4a, except for the change in optics 
from SHAG to LHAG, yielded the results shown by the data points in 
Fig. 4b. The solid line is the prediction of the model based on the 
value of C Q calculated from Eq. 68. The value of the baseline plasma 
ion energy cost was held constant at e* = 57 eV since changing the 
optics should not affect this parameter. The degree of agreement be- 
tween the data points and the curve In Fig. 4b shows clearly tnat the 
model correctly predicts the variation in the plasma ion energy cost 
with the neutral density parameter when the grid system transparency to 
neutrals is changed. 

The same procedure of calculating a now value of the primary 
electron utilization factor from the old value according to Eq. 26 was 
followed for the analysis of the data displayed in Figs. 5a and 5b. 

For the data in Fig. 5a, the thruster was operated with krypton 



Figure 5a. Plasma Ion Energy Cost Curve for Krypton - Configuration I 



Figure 5b. Plasma E-Jjy Cj.t Curve for Small Diameter Ion 
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propellant and SHAG optics at a discharge voltage of 40 v. The primary 
electron utilization factor corresponding to this operating condition 
was calculated using the value of C Q obtained from Fig. 4a together 
with Eq. 26 In the form. 


^o^new 


(O 


o'Kr 


K>Ar 



(C o^F1c. 4a 


(69) 


In this case, both the change In propellant properties and the change 
in discharge voltage must be accounted for. That is, (<Xg) Kr Is the 
total inelastic collision cross section for 40 eV primary electron- 
krypton atom collisions, whereas (o' Q )^ r refers. In this case, to 50 eV 
primary electron-argon atom collisions. The cross section data needed 
in this equation were obtained from de Heer, et. al [48]. The new 
value of the primary electron utilization factor calculated from 
Eq. 69 was C Q = 5.7 (A eq.)" 1 and the corresponding prediction of the 
model is compared to the measured values in Fig. 5a. A new value of 
the baseline plasma ion energy cost (e*) was also required to fit the 
data of Fig. 5a since this parameter Is a function of both the dis- 
charge voltage and the propellant type. This new value was selected 
in the manner suggested previously as the measured plasma Ion energy 
cost at high neutral density levels. These data Indicate that the 
model works as well for operation with krypton propellant as it does 
for argon. 

From Eq. 26, it is seen that the primary electron utilization 
factor depends inversely on the area of the grids through which the 
beam Is extracted, A g . This area may be varied without changing the 
discharge chamber diameter by masking down the screen grid to produce 
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ion beams of different cross sectional areas. In this case, the screen 
grid for the SHAG optics set was masked down from a beam diameter of 
12 cm to oneof 6cm. This four fold reduction in beam area should pro- 
duce a corresponding four fold increase in nrimary electron utilization 
factor. Taking C Q for Fig. 5a and multiplying by four yields the new 
value of C Q = 22.8 (A eq.)" 1 . The model prediction for the plasma ion 
energy cost versus neutral density parameter curve, obtained using this 
value of C Q , is compared to the measured values of these quantities in 
Fig. 5b. Remarkably, the agreement between the model and the experiment 
is excellent. Similar agreement was obtained for operation with argon 
using the masked down grid set. 

Measurements of the plasma ion energy cost variation with the 
neutral density parameter were also made on the thruster configuration 
of Fig. 3b. Some of the results obtained with this configuration are 
shown in Figs. 6a, 6b, 7a and 7b. The data in Fig. 6a was obtained 
with argon propellant, a discharge voltage of 50 v and a range of pro- 
pellant flow rates from 350 to 1500 mA eq. Excellent agreement between 

★ 

the model and the experimental data was obtained by taking e p = 59 eV 
and C Q - 4.0 (A eq.)" 1 . The value of e* = 59 eV, obtained with this 
configuration, agrees well with the value of e* = 57 eV obtained with 
the other configuration for operation with argon at a discharge voltage 
of 50 v. The slightly higher value of ej for the configuration of 
Fig. 3b is believed to be the result of the shorter discharge chamber 
length causing a slightly higher fraction of ion current to go to the 
anod® surface * n d rathode support posts. 
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From the value of C Q in Fig. 6a, a new value of C Q applicable to 
this configuration operating on xenon propellant at a discharge voltage 
of 40 v may be calculated using Eq. 26. The results of this calcula- 
tion, together with the experimental results, are given in Fig. 6b. 

For these data, the neutral density parameter has been corrected for 
presence of doubly charged ions in the beam. Again the value of e* was 
chosen in the manner described above. Clearly, the model agrees well 
with the measured results. Equations 25 and 26 have now been demon- 
strated to describe correctly the effects of variations in the plasma 
ion energy cost resulting from changes in propellant flow rate, pro- 
pellant utilization, effective grid transparency to neutral atoms, beam 
extraction area, discharge voltage, and propellant gas (Ar, Kr and Xe). 
In addition, other experiments indicate that the model correctly pre- 
dicts the change in thruster performance resulting from a change in the 
effective discharge chamber wall temperature [36]. The model has also 
been shown to work wall on line cusp thruster geometries [35]. 

For operation at low discharge voltages, however, the situation is 
quite different as seen in Fig. 7a. The data in this figure were taken 
using argon propellant, a discharge voltage of 30 v and the thruster 
configuration of Fig. 3b. Here, a systematic difference in the data 
taken at different propellant flow rates is observed. Clearly, a single 
equation such as Eq. 25 is not sufficient to explain this behavior if 
C Q and e* are taken to be independent of the propellant flow rate and 
propellant utilization. This systematic difference of the plasma ion 
energy cost curves, measured at different flow rates. Is believed to 
result from changes in the baseline plasma ion energy cost (e*) which 
occur at low discharge voltages. 


These changes, in turn, result from the dependence of the baseline 
plasma ion energy cost on the average energy of a Maxwellian electron 
collected by the anode, (see Eq. 27). The value of e^, however, 
depends on the Maxwellian electron temperature at the anode. Measured 
values of the electron temperature, made using the Langmuir probe po- 
sitioned on the surface of the anode magnet ring, are shown in Fig. 7b. 
These data were taken under operating conditions identical to those 
used for the collection of the data in Fig. 7a. The data in Fig. 7b 
clearly indicate a systematic difference in electron temperature at the 
anode for different propellant flow rates. For a given value of the 
neutral density parameter, these data show that the electron temperature 
at the anode increases with increasing propellant flow rate. Higher 
electron temperatures at the anode correspond to higher values of e* 
because of the increased loss of energy from the plasma in the form of 
Maxwellian electrons. Because the ratio e^/Vg appears in the equation 
for e* , changes in e M become increasingly important at low discharge 
vo 1 tages . 

The systematic difference in the curves of electron temperature 
versus neutral density parameter observed in Fig. 7b is not predicted 
by the model (as given by Eq. 59). Thus, some additional mechanism for 
the transfer of energy from the primary electrons to the Maxwellian 
electrons, othT than that considered in the model, must become im- 
portant at low discharge voltages. This mechanism is believed to be 
the direct thermal Ization of the primary electrons as the result of 
electron-el ectron collisions. The collision cross section for primary- 
Maxwellian electron collisions increases rapidly with decreasing 
primary electron energy. This is consistent with the observation that 
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the separation between the electron temperature curves taken at dif- 
ferent propellant flow rates (such as those in Fig. 7b) Increases with 
decreasing primary electron energy (i.e., discharge voltage). 

Extracted Ion Fraction 

Aside from the plasma ion energy cost* the other parameter which 
has a major affect on thruster performance is the extracted ion frac- 
tion, f D . It is of little use to create ions efficiently in the dis- 
charge chamber if the fraction of these ions extracted into the beam 
is small. The extracted ion fraction, which is defined In Eq. 6 as the 
ratio of the beam current to the total ion current produced, was mea- 
sured for both thruster configurations over the range of operating 
conditions discussed earlier. 

The resulting extracted ion fraction data, obtained on the 
thruster configuration of Fig. 3a, are shown in Figs. 8a and b. In 
Fig. 8a, the value of f g are given versus the neutral density parameter 
for operation with argon propellant at discharge voltages of 30 and 
50 v over a range of propellant flow rates from 500 to 1500 mA eq. 

These data indicate that the extracted ion fraction, while being rela- 
tively independent of neutral density, does tend to be slightly greater 
at lower discharge voltages. Data taken under the same conditions, but 
at a 40 v discharge voltage, fall between the two curves shown in Fig. 8a. 

The extracted ion fractions for argon propellant are compared to 
those for krypton propellant at a discharge voltage of 40 v In Fig. 8b, 
The values of f g for argon are seen to be generally slightly higher 
than those for krypton. 
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The data in Figs. 8a and 8b indicate that the extracted ion 
fraction is relatively independent of the neutral density parameter. 
Other data not presented, indicate that the extracted ion fraction 
sometimes decreases with increasing beam currents for certain thruster 
configurations. This decrease in f g results from a decrease in the 
effective transparency of the screen grid to ions as the plasma density 
is increased. However, the fraction of the total ion current produced 
that is directed toward the grids remains constant. 

Earlier studies [32,33] indicate that the extracted ion fraction, 
or the fraction of ions directed toward the grids, is strongly dependent 
on the magnetic field configuration and thruster geometry. The data in 
this report indicate that fg is also a weak function of the discharge 
voltage and propellant gas, but not a function of the neutral density 
parameter. Consequently, a simple design approach would be to take f„ 
to be a constant for a given discharge chamber design (magnetic field 
configuration), propellant and discharge voltage. Unfortunately, a 
method for calculating the value of the extracted ion fraction from the 
above design data remains to be developed. 

Plasma Properties 

Equation 45 provides a simple expression for the primary electron 
density as a function of the propellant utilization. As mentioned 
earlier, the combination of parameters inside the square brackets in 
this equation may be taken to be roughly a constant for a given pro- 
pellant, discharge roltage and thruster configuration. The value of 
this constant, appropriate for argon propellant, a discharge voltage of 
50 v and the thruster configuration of Fig, 3b, was calculated to be 
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1.0 x 10 10 cm" 3 . Thus, Eq. 45 becomes, In this case, 


n = (1.0 x 10 10 ) (cm -3 ) . (70) 

F n u 

The value of e* , required for this calculation, was taken to be the 
measured value given in Fig. 6a. The value of e* could, however, have 

r 

been calculated using the method described in Appendix A. 

The volume of the ion production region (V-), required by Eq. 45, 
was determined from a computer drawn magnetic field map of the dis- 
charge chamber. This map was created by measuring the magnetic flux 
density and direction at regularly spaced points in the discharge 
chamber. The ion production volume, defined by this map, was taken to 
be the volume in which the magnetic flux density was 0.005 Tesla or 
less. 

The extracted ion fraction (fg) was taken to be its measured value 
obtained in tests with argon propellant, a discharge voltage of 50 v 
and the thruster configuration of Fig. 3b. The transparency of the 
grids to neutral atoms (<t> Q ) was calculated from Eq. 66. Finally, the 
neutral atom velocity (v Q ) was calculated based on an assumed effective 
wall temperature of 400 K. 

Comparison of the calculated values of the primary electron density 
from Eq. 45 (in the form of Eq. 70) with the measured values is given 
in Fig. 9a. The measured primary electron densities were obtained using 
the Langmuir probe positioned along the centerline of the discharge 
chamber. Equations 45 and 70, however, calculate essentially an average 
primary electron desnity. Further, since the centerline primary elec- 
tron density is expected to be higher than the average value, the 
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excellent agreement between the calculated and measured values In 
Fig. 9a is somewhat misleading. This figure clearly indicates, however, 
that Eq. 45 has the correct functional form and that the primary elec- 
tron density is indeed relatively independent of the propellant flow 
rate as predicted by the model . 

A similar comparison between calculated and measured primary 
electron densities is given in Fig. 9b for operation of the same 
thruster configuration, but with xenon propellant and a discharge volt- 
age of 40 v. Again, the same concludsions drawn from Fig. 9a are also 
applicable to thruster operation on xenon as indicated in Fig. 9b. The 
observation that the calculated average values are somewhat greater 
than the measured centerline values probably results from the use of 
the measured value of e* = 44 eV in Eq. 45 rather than the calculated 
value of e* = 36 eV given in Appendix A. 

An expression for the ratio of primary- to -total electron density 
is given by Eq. 48. The combination of parameters in the square brac- 
kets in this equation is approximately a constant for a given propellant, 
discharge voltage and thruster configuration. For operation with argon 
at V D = 50 and the configuration of Fig. 3b, Eq. 48 becomes, 


!!r = 0.027 

n i m(l-n u ) 


(71) 


In making this calculation, the Bohm velocity (v b ) was calculated based 
on an electron temperature of 4 eV. The value of the screen grid 
transparency to ions ($.) was determined experimentally by measuring 
the ion current to the screen grid and to the beam. The transparency 
is then the ratio of the beam current to the sum of the screen grid and 
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beam currents. The measured screen grid transparency to Ions was 
approximately 0.8 compared to the physical open area fraction of the 
screen grid which was approximately 0.68. 

Comparison of Eq. 48 (in the form of Eq. 71) with the experimental- 
ly measured values is given in Fig. 10a. This figure indicates that for 
a given discharge voltage, propellant and thruster configuration the 
primary- to- total electron density ratio is only a function of the neu- 
tral density parameter. In addition, it indicates that Eq. 71 has the 
correct functional form for the variation of this ratio with the 
neutral density parameter. 

For operation of the same thruster configuration with xenon at 
V Q = 40 v, the results shown in Fig. 10b were obtained. The solid line 
in this figure corresponds to Eq. 48, where the constant was calculated 
using values of the parameters in Eq. 48 that are appropriate for xenon 
propellant and a discharge voltage of 40 v. Figures 10a and b indicate 
that Eq. 48 correctly accounts for changes in the propellant gas, the 
discharge voltage, the neutral density parameter and the propellant 
flow rate. 

Comparison of calculated and measured Maxwellian electron tempera- 
tures is given in Fig. 11. The calculated electron temperatures were 
obtained using Eq. 59. For operation with xenon at Vq = 40 and the 
thruster configuration of Fig. 3b, Eq. 59 becomes. 


q + = 7.64 x 10~ 14 

0 226 m(l-n u )-l 


(m 3 /s) 


(72) 


where the appropriate cross section data were obtained from References 
49 and 50. Equation 72 gives the value of the Maxwellian ionization 
rate factor for the production of single ions from neutral atoms as a 
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function of the neutral density parameter. The rate factor Q* Is 
function of the Maxwellian electron temperature, and Is given by. 





a 


(73) 


Using Eg. 73, the rate factor Q* may be plotted as a function of the 
electron temperature as shown In Fig. 12. 

To determine the electron temperature as a function of the neutral 
density parameter, the following procedure is used. First, the value 
of the rate factor Q* is calculated for a given value of m(l-n u ) using 
Eq. 72. This value of Q* is then used to enter the curve in Fig. 12 
from which the corresponding electron temperature is determined. Re- 
peating this procedure generates the curve of Maxwellian electron 
temperature versus neutral density parameter shown as the solid line 

m Fig. 11. The agreement between the calculated and measured electron 
temperatures is considered to be quite good. 

In addition to the Langmuir probe positioned on the thruster 
centerline, a second probe was positioned on the upstream magnet ring. 
This probe was used to measure the temperature of the Maxwellian elec- 
trons and the energy of the primary electrons reaching the anode. 

Because the probe was positioned in a region of very high magnetic flux 
density, the electron temperatures determined from the probe traces 
obtained with this probe are probably only equal to the electron tem- 
perature resulting from motion along the magnetic field lines. The 
electron temperature corresponding to motion normal to the field lines 
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may or may not be the same. For simplicity, the^e temperatures are 
assumed, in this work, to be the same. 

The measured electron temperatures at the anode were found to be 
less than the measured centerline temperatures by approximately a 
factor of 2/3. The correlation of electron temperature at the anode 
with the centerline electron temperature is given in Fig. 13. The data 
in this figure were obtained with the thruster configuration of Fig. 3b 
using both argon and xenon propellants. The solid line in this figure 
has a slope of 2/3. This observation that the electron temperature at 
the anode is approximately 2/3 of the centerline temperature is used in 

Appendix A for the calculation of e* . 

P 

The probe traces obtained with the probe positioned on the magnet 
ring surface also indicated the presence of primary electrons. For 
operation at a discharge voltage of 50 v, these primary electrons had 
an energy of approximately 50 eV as the model suggests they should for 
this case where a refractory cathode is used (V c = 0). Similarly, for 
operation at V Q = 40 v, the measured primary electron energy at the 
anode surface was ^ 40 eV. These probe traces provide direct evidence 
of the loss of primary electrons through the magnetic field cusp to the 
anode. Further, the loss rate of these primary electrons was seen to 
Increase as the neutral density parameter decreases as the model 
predicts. 

A comparison between the calculated (from Eq. 66) and measured 
values of the doubly-to-slngly charged ion beam current ratio Is given 
in Fig. 14. The measured values are given by the open symbols and the 
calculated values by the solid symbols. The propellant utilization 
efficiencies in this figure have been corrected for the presence of 











doubly charged ions according to Eq. B-5 given in Appendix B. The 
ratios Jg + / Jg in Fig. 14 were taken to be equal to the values of 
doubly-to-singly charged ion current density measured on the thruster 
centerline. The values calculated from Eq. 66, however, correspond to 
values of this ratio averaged over the entire grid area. Since the 
centerline ratio is expected to be greater than the average value, the 
data points in Fig. 14 are believed to be shifted up and to the left 
relative to the locations of the corresponding average values. The 
amount of this shift is unknown. This shift, however, is believed to 
account for at least part of the difference between the measured and 
calculated values. 

In calculating the values of Jg + /Jg usin 9 E< 1* the collision 
cross section data required for the calculation of the parameters Q* 
and Pg were taken from Reference 50. The values of the parameters 
Q 0 ++ . 3+*^ P o + > p + + were obtai 'ned from Reference 11, in which the colli- 
sion cross sections were calculated using the classical cross section 

model of Gryzinski [51]. The use of Gryzinski cross section data in 
the calculations of Jg + /Jg might also account for some of the difference 
between the measured and calculated values. The values of the Max- 
wellian electron temperature and primary- to-Maxwel 1 ian electron density 
ratio, required as inputs to Eq. 66, were calculated using Eqs. 48, 49 
and 59. 

Finally, it is noted that, the values of Jg + / Jg calculated from 
Eq. 66 are very sensitive to the value of the electron temperature used 
in the calculations. Equation 66 does, however, indicate the correct 
functional dependence of the doubly-to-singly charged beam current 
ratio on thruster operating conditions. 


V. MODEL APPLICATIONS 


Thruster Design 

Equation 28 (repeated below) provides a single equation describing 
the performance curve of a given thruster design 


£n - 


f C V D 




This equation, along with the data presented in Chapter IV, suggest that 
the performance of any thruster design depends on the values of four 
discharge chamber configuration controlled parameters: e*, C Q , fg and 
f c ; and two operating parameters, m and Vp. The effect of these 
parameters on performance was investigated analytically through Eq. 28 
by choosing the following set of values as the standard set, then vary- 
ing them one at a time to determine their effect on the traditional 
performance curves. 


Table 1 

Standard Configuration Parameters 
m = 1000 mA eq. 
f B =0.6 
f c =0.1 
V L , = 50 v 
C Q = 3.0 A eq. 
e* = 50 eV 
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Figure 15a shows the effect of the extracted ion fraction on 
performance. As expected, this parameter has a strong effect on the 
performance. Changes in fg shift the performance curve up or down, but 
do not substantially change its shape. Clearly, it is desirable to 
have fg be as large (near unity) as possible. 

The effect of the primary electron utilization factor (C ) on 
performance is given in Fig. 15b. This parameter also has a strong 
effect on the performance. Indeed, it is this parameter which primarily 
determine the shape of the performance curve, with larger values of C Q 
corresponding to improved performance and curves with more sharply 
defined "knees." The definition of C Q given in Eq. 26 suggests a num- 
ber of ways in which the value of C Q may be increased. For example, it 
may be increased by using a propellant gas characterized by a larger 
inelastic collision cross section a^, and a larger atomic mass (result- 
ing in a lower neutral velocity v Q ). The parameter C Q may also increase 
by decreasing the grid transparency to neutrals, $ . This must be done, 
of course, without increasing the accelerator impingement current. For 
thruster designs with non-uniform beam profiles, Tailoring the accele- 
rator grid hole size to match the radial current density profile might 
be a useful way to minimize <t> Q . Also, three grid systems might be ex- 
pected to have smaller values of 4 > Q than two grid systems. 

Most importantly, C Q may be increased by increasing the primary 
electron containment length {Z ). This length corresponds to the 

c 

average distance a primary electron would travel in the discharge 
chamber before being collected by an anode surface provided it had no 
inelastic collisions. As mentioned earlier, the primary function of the 
magnetic field in all thruster designs is to increase this length. In 






cusped field thrusters, primary electrons are lost to the anode through 
the cusps. Thus, z may be increased by decreasing the number of cusps 
at anode potential or increasing the flux density at the cusp, but only 
up to a point. Reductions in effective anode cusp area below a certain 
limit will result in unstable operation of the discharge [22]. 

Equation 26 also suggests that the primary electron utilization 
factor may be increased by masking down the area of the grids through 
which the beam is extracted, A g . However, decreasing A g in this manner 
will lead to a large reduction in the extracted ion fraction, and, 
therefore, an overall reduction in performance. 

The effect of propellant flow rate on performance is shown in 
Figs. 16a and b. In general, higher flow rates produce better perform- 
ance. The maximum flow rate at which the thruster can be operated, 
however, is limited by the ability of the accelerator system to extract 
the ion current directed toward it. The effect of flow rate on perform- 
ance is less dramatic for thruster designs characterized by larger 
values of C Q as shown in Fig. 16b. High values of C Q should, therefore, 
be particularly desirable in thrusters designed to be throttled. 

The effect of ej is merely to shift the performance curves up or 
down. The amount of the shift increases for smaller values of f g . For 
thrusters that use hollow cathodes, the efficiency of the cathode opera- 
tion (characterized by V c in Eq. 27) has a strong effect on the value 

of e* . In addition, higher values of V D , in general, produce smaller 

P 

values of e* . High discharge voltages, however, are undersirable from 
thruster lifetime considerations, Consequently, a trade off between 
thruster performance and thruster lifetime is necessary. 
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The effect of cathode operation on thruster performance can be 
Investigated analytically using the performance model. Hollow cathode 
operation is reflected in the value of the parameter V c . The lower 
the value of this parameter the more efficiently the cathode is oper- 
ating. At a fixed discharge voltage, increasing the value of V^ 
degrades the thruster performance in two ways. First, it increases 
e* directly through its appearance in Eq. 27, and second, it decreases 
the energy of the primary electrons, which increases the value of the 
collisional loss parameter e Q and decreases the value of C Q . All of 
these affects are accounted for when e* is calculated according to the 
procedure described in Appendix A and the primary electron energy is 
given by V Q -V C . 

The effect of V^. on the baseline plasma ion energy cost was calcu- 
lated for xenon propellant using the method described in Appendix A and 
assuming a discharge voltage of 30 v. The results of these calculations 
are given in Table 2. 


Table 2 

Effect of V c on e* 


( Vol ts ) 

(Volts) 

(Volts) 

(eV) 

(A eq.)* 1 

30 

0 

30 

36 

8.0 

30 

2.5 

27.5 

42 

7.9 

30 

5.0 

25.0 

48 

7.6 

30 

7.5 

22.5 

55 

7.2 

30 

10.0 

20.0 

66 

6.5 

30 

12.5 

17.5 

79 

5.2 
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The values of tp and C Q , from Table 2, were then used In the equation 
for thruster performance (Eq. 28) along with the values: m = 1000 mAeq., 
fg = 0.5, f^ = 0.1, and Vg = 30 v. The performance curves generated in 
this manner are shown in Fig. 17. This figure indicates that, at a dis- 
charge voltage of 30 v, relatively small changes in can produce 
substantial changes in thruster performance. For example, Increasing 
V c from 10 to 12.5 v causes an increase of approximately 30 eV/beam ion 
in the performance curve. 

Work done by Siegfried [52] on inert gas hollow cathodes indicated 
that a 0.6 v Increase in the surface work function of the cathode in- 
sert could produce an increase in internal cathode plasma potential 
from 8 to 12 v. The value of would be expected to Increase by this 
same amount. In addition, Siegfried notes that there is apparently a 
greater sensitivity of the insert to depletion or contamination of the 
low work function material for opjration with argon or xenon as compared 
to operation with mercury. Consequently, on the basis of these consid- 
erations, one would expect the performance of inert gas ion thrusters 
operating at low discharge voltages to be quite sensitive to the con- 
dition and operation of the hollow cathode. 

Changes in the fraction of ion current going to cathode potential 
surfaces (f c ) also tend to shift the performance curves up or down. 

Recent trends [23,24] in thruster design to operate the discharge 
chamber body at anode rather than cathode potential serve to reduce the 
value of f c , and consequently improve the performance. In addition, 
elimination of the separate cathode discharge region in these designs 
serves to decrease both f^ and V^, which again improves the performance. 
It should be noted, however, that an increase in the extracted ion 
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figure 17. Effect of Cathode Operation on Performance for Xenon 
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fraction fg Is also believed to be an Important factor In the perform- 
ance Improvement observed In these new thruster designs. 


Thruster Scalint 


The discharge chamber model also suggests performance changes that 


might be expected from the scaling of thruster designs to different 
discharge chamber diameters. This may be seen by examining the product 
of the primary electron utilization factor and propellant flow rate. 


This product is a dimensionless quantity and is given the symbol B Q . 



4 Ve m 

e Vo A g 


(74) 


For constant, average beam current densities, thruster scaling should 
be accomplished such that the ratio of mass f*ow rate to active grid 
area, m/A , is constant. Thus, for the same propellant and grid trans- 
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parency to neutrals, Eq. 74 suggests that the discharge chamber perform 
ance depends only on the primary electron containment length, It i 
believed that this length should increase in approximately direct pro- 
portion to thruster diameter, and this suggests that larger diameter 
thrusters should be more efficient than smaller ones. This is not the 


whole story, however. The effect of thruster diameter on the extracted 
ion fraction must also be considered. The above conclusion remains 


true for thruster scaling accomplished in such a way that fg remains 
constant. 


Neutral Loss Rate 

Rewriting Eq. 22 yields an expression for the neutral atom loss 
rate from the discharge chamber through the g*"1ds. 




, je* 




78 


n Q = m - J B = »0-n u ) * 71 

Substituting this Into Eq. 28. and solving for the neutral loss rate 
yields, 


This equation gives the value of the neutral atom loss rate as a func- 
tion of the beam ion energy cost. For a specified thruster geometry 
and discharge voltage, the design parameters C 0> e*. f B . f c and V„ may 
be taken to be approximately constant. Thus, since the propellant mass 
flow rate doesn't appear on the right-hand-side of Eq. 76, this equa- 
tion predicts that the neutral loss rate n 0 Is Independent of the flow 
rate at a constant value of the energy cost per beam Ion U B ). This 
same conclusion was reached originally by Kaufman [53] In his constant 
neutral loss rate theory developed for low magnetic field strength dls- 
charge chamber designs. 


i - f - 5 r r 
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Thruster Testing without Be am Extraction 
It Is often desirable to compare the performance of different dis- 
charge chamber designs that have been operated without Ion beam extrac- 
tion. Operation without beam extraction Is often more convenient and 
generally requires the use of smaller vacuum test facilities than 
operation with beam extraction. However, data obtained under these 
conditions must be Interpreted carefully. It has been observed that 
the performance (eV/beam Ion) extrapolated from thruster operation 
without beam extraction Is generally significantly better than the 
performance measured with beam extraction [24,54.55]. This difference 
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in performance may be understood In context of the performance model 
developed in this report. 

Thruster operation without beam extraction should be characterized 
by higher discharge chamber neutral densities than operation at the 
same propellant flow rate with beam extraction. This is because, 
without beam extraction most of the propellant leaves the discharge 
chamber in the form of neutral atoms at the neutral atom thermal 
velocity. With beam extraction, however, most of the propellant leaves 
in the form of ions at the Bohm velocity. Therefore, since the total 
propellant flow rate leaving the thruster must be the same in each case, 
the loss rate of neutrals is smaller with beam extraction than without 
it, implying lower neutral densities. This is primarily a consequence 
of the changing effective transparency of the screen grid. Without 
beam extraction, the effective screen grid transparency is very small 
[56] since ions tend to be focused onto the grid webbing. With beam 
extraction, the ions tend to be focused away from the screen grid 
webbing and through the grid apertures. In any case, operation at 
higher neutral densities, for the same thruster configuration, trans- 
lates into lower plasma ion energy costs according to Eq. 25. 

The performance model developed herein may be used to make more 
meaningful extrapolation of data taken without beam extraction to opera- 
tion with beam extraction. To do this, it is necessar> to experiment 
ally determine the values of the parameters e*. C Q , f B and f c . The 
parameters e* and C Q are roughly independent of the neutral density, 
thus, shouldn't change depending on whether or not a beam is extracted. 
For hollow cathode equipped thruster designs, however, the baseline 
plasma ion energy cost includes the parameter V c . which is an indication 


* 
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of cathode performance. The value of V c may be a function of the dis- 
charge chamber neutral density and consequently could cause the value 

of e* to change for operation with and without beam extraction. The 
P 

magnitude of the change in Vq, if any, is unknown. 

To determine the values of e* and C Q for operation without beam 
extraction, some method of estimating the total ion current produced as 
a function of the discharge chamber neutral density must be used. This 
is probably most easily accomplished using several ion current probes 
positioned at the walls of the discharge chamber in the manner suggested 
by Poeschel [54]. Ideally, these data should be acquired over a range 
of neutral densities at a constant discharge voltage. 

Once the total ion current produced is known as a function of the 
neutral density, a curve similar to fig. 4a may be generated, where the 
values of ep are calculated using Eq. 2. From these data, the values 
of £* and C Q may then be determined as those which give the best fit 
of Eq. 25 to the data. Note, e* may also be calculated according to 

r 

the procedure outlined in Appendix A. 

Finally, the values of fg and fQ may be determined using the above 
probe data and guessing a value of the effective screen grid transpar- 
ency appropriate for the beam extraction condition. Once the values of 
e*, CL, f D and f r have been determined, the performance curve for the 
beam extraction condition may be approximated using Eq. 28, for any 
desired total propellant flow rate. 
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Space Propulsion Mission Analysis 
The performance model, developed in this report, is particularly 
useful for examining the impact of thruster performance on space pro- 
pulsion missions. It should also be useful for studies making compari- 
sons of the capabilities of different propulsion systems such as those 
done in References 57 and 58. The following analysis is intended as an 
illustration of the kinds of things that can be done using this thruster 
performance model. Consequently, the orbit transfer mission considered 
here has been greatly simplified. Solving a more accurate and compli- 
cated orbit transfer problem would not provide any additional insight 
into the applicability of the model to mission analysis, and may even 
tend to conceal the desired illustration. In the analysis that follows, 
the effect of thruster performance on the maximum payload fraction ob- 
tainable for a low earth orbit to geosynchronous earth orbit transfer 
mission, with a characteristic velocity of 6000m/ s, is investigated. 

The rocket equation (Eq. 77) gives the ratio of final spacecraft 
mass (M f ) to initial mass (M i ) as a function of the characteristic 
velocity for the mission (aV), the thruster exhaust velocity (u) and 
the thruster propellant utilization efficiency (n^) » 



-AV/n u u 


(77) 


Assuming, for simplicity, that the final mass consists only of the 
payload mass (M £ ) and the mass of the power plant (M g ) then Eq. 77 can 


o 


be written as, 



The mass of the power plant is proportional to the power generated. 


Mg = a P g , (79) 

where a is the specific power plant mass (Kg/W), and P g is the 
generator power (W). The generator power may be given as. 
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where is the thruster electrical efficiency and m p is the total 
propulsion system propellant mass flow rate (kg/s). For a constant 
exhaust velocity and propellant efficiency, Eq. 80 may be integrated to 
obtain, 
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where M p is the initial propellant mass and t is the total mission 
time. Combining Eqs. 78-80 and recognizing that, 
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Considering only discharge chamber losses, the thruster electrical 
efficiency may be approximated by, 


= 1 
1+e B /V N 


(84) 


where V N is the net accelerating voltage. The net accelerating voltage 
is related to the exhaust velocity by, 
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Combining Eqs. 83-85 yields. 
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Equation 86 gives the payload mass fraction as a function of: the 
characteristic velocity, the exhaust velocity, the mission time, the 
power plant specific mass, the charge- to-mass ratio of the beam ions, 
the average beam ion energy cost and the propellant utilization effici 
ency. The average beam ion energy cost is related to the propellant 
utilization efficiency through the equation developed in the thruster 
performance model , 


e B = 


f B n-e 




Vo 


(87) 


where B Q is a dimensionless quantity describing the utilization of 
primary electrons in the discharge chamber and is given by Eq. 74. 

For a given characteristic velocity, mission time, power plant 
specific mass, and propellant utilization efficiency, the optimum 
exhaust velocity corresponding to the maximum payload fraction, for 
these conditions, can be determined from Eq. 86. This is accomplished 
by setting the derivative of Eq. 86 (with respect to u) equal to zero 
and numerically solving the resulting equation for u. Substituting 
this value of the exhaust velocity back Into Eq. 86 gives the value of 
the optimum payload fraction under the specified conditions. This pro 
cedure may be repeated for different values of propellant utilization. 
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where, for each value of n u » the corresponding beam Ion energy cost is 

calculated from Eq. 87. In this manner, a curve of optimized payload 

fraction versus propellant utilization efficiency may be generated for 

specified values of the characteristic velocity, mission time, power 

plant specific mass and thruster performance parameters, e* f D , f., 

P B C 

B q and Vg. 

An example of this is shown as the solid line in Fig. 18 for a 
characteristic velocity of 6000 m/s, a mission time of 200 days and a 
specific power plant mass of 40 kg/kW. In addition, the thruster per- 
formance parameters were taken to be e* * 50 eV, f R ■ 0.5, f r = 0.1, 

B 0 * 5 and Vg * 40 v with argon as the propellant. Figure 18 Indicates 
that the curve of optimized payload fraction versus propellant utiliza- 
tion efficiency goes through a maximum. The propellant utilization 
efficiency corresponding to this doubly maximized payload fraction In- 
dicates the location on the thruster performance curve at which the 
thruster should be operated In order to truly maximize the payload 
fraction. In summary, the doubly maximized payload fraction was de- 
termined by optimizing the exhaust velocity to produce the optimum pay- 
load fraction for a given propellant utilization efficiency and 
subsequently selecting the propellant utilization to obtain the maximum 
payload fraction from this set of optimized payload fractions. 

To Illustrate more clearly what is going on, the Initial propellant 
mass fraction ( M p Q /M^) and power plant mass fraction (M g /M^) are plotted 
along with the optimized payload fraction in Fig. 18. Note that the 
sum of the three curves in this figure, at any propellant utilization, 
is equal to unity. As the propellant utilization increases the initial 
propellant mass fraction decreases as expected. In addition, the power 
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plant mass decreases Initially, due to a decreasing optimum exhaust 
velocity, then Increases dramatically at high propellant utilization 
efficiencies, due to the rapid increase in the average beam ion energy 
cost. For the set of conditions chosen for the curves in Fig. 18, the 
doubly maximized payload fraction occurs at a propellant utilization 
of 0.73. 

The effect of the specific power plant mass on the optimized pay- 
load fraction vs. propellant utilization efficiency curve is shown in 
Fig. 19a. It is interesting to note that, at very low values of a the 
optimized payload fraction is relatively insensitive to large changes 
in the propellant utilization. In addition, the propellant utilization 
at which the doubly maximized payload fraction occurs is a function of 
a. That is, the point at which one should operate on the performance 
curve of a given thruster depends on the specific power plant mass for 
the spacecraft. This is illustrated in Fig. 19b, where the optimum 
propellant utilization efficiences for the specified mission parameters 
are indicated on the performance curve used in generating the curves in 
Fig. 19c. It is also Interesting to note that, the optimum propellant 
utilization for the o 8 1 kg/kW case occurs well past what might ordin- 
arily be considered the "knee" of the performance curve. In a similar 
manner, the mission time can be shown to effect the optimum propellant 
utilization, with longer trip times yielding higher optimum propellant 
utilizations. 

So far, only one thruster performance curve has been considered and 
it has been observed that the optimum propellant utilization efficiency 
corresponding to the doubly maximized payload fraction Is a function of 
the specific power plant mass and the mission time. The effect of the 
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thruster performance parameters on the optimum payload vs. propellant 
utilization curves will now be Investigated. 

Figure 20a indicates the effect of the parameter B Q on payload 
fraction curves. This parameter primarily determines the shape of the 
thruster performance curve. As expected, larger values of B Q yield 
higher maximum payload fractions. The optimum propellant utilization 
efficiencies are indicated, for each performance curve in Fig. 20 b, by 
the vertical lines that intersect the curves. 

The effect of the extracted ion fraction on the payload fraction 
curves is shown in Fig. 21a. Not surprisingly, higher extracted ion 
fractions produce larger maximum payload fractions. Again, the optimum 
propellant utilizations are indicated for each performance curve in 
Fig. 21b. 

Finally, the effect of the propellant gas is indicated in Fig. 22a. 
To generate the curve labelled "argon" in this figure, the following 
thruster performance parameters were used: f g = 0.5, fp = 0.1, Vp * 

40 v, B q = 5.0 and e* = 50 eV. For krypton and xenon propellants, the 
values of f g and f^ were held constant, but, the values of Vp. B Q and 
e* were changed to 30 v, 8.2 and 49 eV for krypton, and 30 v, 14.6 and 
42 eV for xenon. The changes in the values of B Q were calculated using 
the equations developed in the thruster performance model. The changes 
in the values of e* are consistent with measured changes in this param- 
eter for a given thruster design operated on the three propellant gases. 
Clearly, xenon propellant is superior to either argon or krypton pro- 
pellants. The performance curves, along with the corresponding optimum 
propellant efficiencies, are given in Fig. 22b for the cases where 
these three propellants are used. 
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Figure 22b. Effect of Propellant on the Performance Curve and the 
Optimum Propellant Utilization 


VI. CONCLUSIONS 


A simple thruster performance model, that has led to an improved 
understanding of ion thruster discharge chamber operation, has been 
developed. This model describes ion thruster performance in terms of 
four parameters: the plasma ion energy cost (e p ), the extracted ion 
fraction (fg), the ion current fraction to cathode potential surfaces 
(f c ) and the discharge voltage ( Vp) . 

The equation developed to describe the behavior of the plasma ion 
energy cost agrees with the results of a variety of experiments. This 
equation provides an expression for the functional dependence of the 
plasma ion energy cost on the neutral density parameter, m(l-n u ). 
Experiments indicate good agreement between the predicted functional 
form of the model and the experimental data. These experiments also 
suggest that the primary electron utilization factor (C Q ) and the base- 
line plasma ion energy cost (ej) are independent of the neutral density 
parameter under many conditions. The model correctly predicts the 
variation In plasma Ion energy cost for changes in: propellant gas, 
grid transparency to neutral atoms, beam extraction area, discharge 
voltage and effective discharge chamber wall temperature. The model is 
applicable to both ring and line cusp thruster designs. 

Measurements of the extracted Ion fraction indicate that this 
parameter Is relatively Independent of the neutral density parameter. 
The extracted ion fraction does, however, appear to be a function of 
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the thruster geometry and magnetic field configuration, and to a lesser 
extent the propellant gas and the discharge voltage. No effective 
method for calculating the extracted Ion fraction based on these con- 
siderations Is available at the present time. 

The values of several discharge chamber plasma properties can be 
easily calculated, as a function of thruster operating conditions, 
using the equations developed in the model. These properties include: 
the primary electron density, the primary- to-total and primary- to- 
Maxwellian electron density ratios and the Maxwellian electron tempera- 
ture. The values of these properties, calculated by the model, agree 
well with the corresponding experimental data. An equation was also 
developed to calculate the doubly-to-singly charged ion beam current 
ratio. Experiments indicate only fair agreement between this equation 
and the measured data. This is believed to be partly the result of the 
technique used to measure the doubly-to-singly charged ion beam current 
ratio, and partly the fact that the calculated values of this ratio are 
very sensitive to the input electron temperature. 

For design purposes, thruster performance may be taken to depend 
on only four physical parameters (the primary electron utilization 
factor C Q , the baseline plasma ion energy cost e*, the extracted ion 
fraction f g and the cathode potential surface ion fraction f £ ) and two 
operating parameters (the propellant flow rate m and the discharge 
voltage V g ). Improved performance should be characterized by large 
extracted ion fractions, long primary electron containment lengths 
(t e ), small effective grid transparencies to neutral atoms U 0 ) and 
operation at high propellant flow rates. The loss of primary electrons 
to the anode determines the shape of the performance curve to first 
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order. The model suggests that thruster designs characterized by 
large values of C Q should allow efficient throttling. In addition* it 
suggests that hollow cathode efficiency becomes increasingly Important 
to the discharge chamber performance at low discharge voltages. 

The thruster performance model can be very useful for mission 
analysis calculations. The model allows one to easily identify the 
optimum propellant utilization at which the thruster should be operated 
for a given mission. In addition, the model makes it easy to determine 
which changes in thruster design or operating parameters have the 
greatest effect on the payload fraction and/or mission duration. 

Finally, a simple direct method for the calculation of the base- 
line plasma ion energy cost has been developed. Calculations of Ep 
agree well with the measured values. 

Suggestions for Future Work 

Two of the parameters identified in this investigation as having 
strong effects on thruster performance were the primary electron con- 
tainment length and the extracted ion fraction. At the present time, 
however, these parameters cannot be calculated directly given only the 
discharge chamber/magnetic field design, propellant gas and discharge 
voltage. Consequently, future work should focus on the development of 
a method by which the primary electron containment length and the ex- 
racted ion fraction can be calculated based only on this Information. 

The experiments described in this investigation were all performed 
using a thruster equipped with a filament cathode. The model suggests, 
however, that the performance of a thruster equipped with a hollow 
cathode can be strongly Influenced by the operation of the hollow 
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cathode. Consequently, there Is a need to verify the predictions of 
the model on a hollow cathode equipped Ion thruster. 

Finally, It Is of Interest to test the validity of the model on 
the older style, low magnetic field strength thruster designs, such as 
the SERT II [63] or the J-Serles [64] thrusters. These thruster designs 
differ from the cusped magnetic field designs considered In this Inves- 
tigation in that electrons. In the discharge chamber plasma, can only 
reach the anode surface by crossing magnetic field lines. In the 
cusped field designs, however, the majority of electrons are believed to 
reach the anode by going along the field lines. The mechanism for elec- 
tron loss in the low magnetic field thruster designs is sufficiently 
different from that in the cusped field designs that experimental veri- 
fication of the model on these low field strength designs is required. 
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APPENDIX A 


THEORETICAL CALCULATION OF THE BASELINE 

PLASMA ION ENERGY COST, e* 

P 

The definition of the baseline plasma Ion energy cost (e*) Is given 

r 

by Eq. 27, and Is repeated here for convenience. 
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The parameter e Q In Eq. A-l accounts for the energy that Is expended In 
Ionization and excitation reactions and Is defined In Eq. 15 as. 
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The brackets In this equation represent the enclosed product averaged 
over the entire electron speed distribution function, l.e.. 
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where F(v g ) represents the entire electron distribution function. 

For a plasma with an electron population characterized by a 
Maxwellian distribution of temperature T^, and a monoenergetlc 
(primary) group of energy Ep, Eq. A-2 becomes, 


(A-4) 
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where < >m represents the enclosed product averaged over the 
Maxwellian speed distribution function. 

The term under the summation sign in Eq. A-4 may be approximated 
by considering only a single equivalent lumped excited state character- 
ized by a total excitat 4 i collision cross section a and a lumped 
excitation energy U ex . For rare gases, U 0X may be approximated by [38], 


u e* ' J < u ,i + • 
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where U £ Is the lowest excitation energy level. Using this lumped 
excitation approximation Eq. A-4 becomes. 
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Substituting Eq. A-6 Into A-l yields the following expression for e* 
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The value of e* may be easily calculated using Eq. A-7 for a given 
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Maxwellian electron temperature, prlmary-to-Maxwelllan electron density 
ratio, discharge voltage and value of V £ . In addition, the average 
energy of a Maxwellian electron lost to the anode (e M ) must be known. 
Reference 19 gives this energy as, 

e M s2T A + V A ‘ ( A - 8 ) 

where Is the electron temperature at the anode and is the 
difference between plasma potential and anode potential. Experiments 
presented earlier Indicated the electron temperature at the anode was 
approximately 2/3 of the centerline temperature. Therefore, e M will 
be taken as, 

e M ’ 2 <! V + \ • (A-9) 

Langmuir probe measurements also Indicate that for the thruster config- 
urations tested In this Investigation was always approximately 2v, 
thus, this value of will be used in these calculations together with 
Eq. A-9. 

For xenon, the values of the total excitation collision cross 
section required by Eq. A-7 were taken from Reference 49. Ionization 
cross section data were obtained from Reference 50. These data, along 
with polynomial curve fits used to facilitate interpolation between 
data points are given in Figs. A-la and A-lb. 

With these data, Eq. A-7 was used to calculate the values of e* as 
a function of electron temperature with the primary-to-Maxwelllan 
electron density ratio as a parameter. The results of these calcula- 
tions are given in Fig. A-2 where the values Vq=40v, V c =0v and V^=2v 
were used. From this figure, it is seen that the baseline plasma ion 
energy cost can vary over a wide range of values, and it is not clear 
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Fiaure A-la. Ionization Collision Cross Section for Xenon (Ref. 50 ) 



Fiqure A-lb. Total Excitation Collision Cross Section for Xenon 
(Ref. 49) 
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at this point which value or values of e* are appropriate. This diffi- 
culty arises because the Maxwellian electron temperature and prlmary-to- 
Maxwelllan electron density ratio may not be selected Independently as 
was done for the calculations In Fig. A-2. 


To resolve this difficulty a second equation for e* is required. 
This second equation for e p may be derived in the following manner. 
Equation 40 provides an expression for the ratio of Ion current pro- 
duced by primary electrons to the total ion current produced, i.e.. 


i_ . 

j p 
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The total Ion current produced, however. Is given by Eq. 52 as, 


J p ■ "o n p e V+ * * n oV <0+, e > H ¥ • (A -”> 

and the ion current produced by primary electrons is given by Eq. 30 as, 

j' = n n ev n (i, V- (A-12) 

p o p p + 

dividing Eq. A-12 by A-11 yields. 
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Equating Eq. A-13 to A-10 and solving for e* yields the desired second 
equation for the baseline plasma ion energy cost. 
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The appropriate values of e* , determined from the correct corres- 
ponding electron temperatures and density ratios, may be found by 
solving Eqs. A-7 and A-14 simultaneously for ej and the electron 
temperature for specified values of the primary- to-Maxwell Ian electron 
density ratio. This procedure is most easily accomplished graphically, 
as shown in Fig. A-3, where the Intersection of the curves corresponding 
to the same values of rip/r^ gives both e* and the electron temperature. 
The locus of these intersection points Indicates the variation of the 
baseline plasma ion energy cost with the electron temperature and 
primary-to-Maxwell ian electron density ratio. Figure A-3 Indicates 
that, under the assumptions used for these calculations, e£ does not 
vary substantially over wide variations in. electron temperature and 
n p / ^ r Vl ' ^ 1S a 9 rees with the experimental observation that e* is a 
constant for operation with xenon propellant at a discharge voltage of 
40 v. 

The measured value of e* is shown along with the calculated values 
in Fig. A-3. The experimental value, at low electron temperatures, is 
seen to be generally slightly higher than the calculated values. The 
most likely explanation for this is a systematic measurement error of 
the total ion current produced as discussed in Appendix B. Also Indi- 
cated in Fig. A-3 is the theoretical maximum value of e* for this case 
(as calculated from Eq. 41). Finally, it should be noted that the 
above results are somewhat sensitive to the choice of e u which is 
assumed here to be given by Eq. A-9. 
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APPENDIX B 
ERROR ANALYSIS 


The experiments described in this investigation require the 
measurement of several quantities. These quantities may be separated 
into four groups according to the technique required to make the mea- 
surement. These four groups are: 

1. Measurement of e„, e? and f D . 

P P o 

2. Measurement of the neutral density parameter, m(l-n u ) 

3. Measurement of the plasma properties: n^ , n p , and T^ . 

4. Measurement of doubly charged ions. 

The errors associated with the measurements of the quantities in each 
of these groups will be discussed separately. In each case, both 
systematic and random errors will be discussed. 

Measurement of e p , e* and f B 

Determination of the quantities in this group requires the measure 
ment of the thruster electrical operating parameters including: the 
discharge voltage, discharge current and beam current. In addition, 
the total ion current produced (J p ) must be measured. The total ion 
current produced is given as the sum of the Ion currents leaving the 
plasma. I.e., 

J p = J B + J C + J A + J 1mp ‘ 

The accelerator grid Impingement current is generally less than 


no 


one percent of the beam current, and can probably be neglected. In this 
Investigation, however, the impingement current was always included in 
the total ion production current for completeness. 

Measurements of the thruster electrical operating parameters were 
made using the ion source instrumentation shown schematically in 
Fig. B-l . The measurement of the ion current to cathode potential sur- 
faces (Jq) was accomplished by biasing these surfaces 30v negative of 
cathode potential. With this bias, electron collection at these surfaces 
was eliminated allowing the incoming ion current to be measured. The 
ion current to anode potential surfaces (J A ) and to the cathode support 
posts could not be measured with this technique. The measured ion pro- 
duction current, consequently, does not include these currents. 

Three potentially significant systematic errors associated with 
the measurements of J p (and thus e p and f g ) have been identified. These 
are: 

1. Neglecting the ion current to anode potential surfaces and the 
cathode support posts. 

2. Secondary electrons emitted by ions striking the negatively 
biased surfaces. 

3. The presence of doubly charged ions. 

The most serious of these is the neglect of the ions reaching the 

anode surfaces and the cathode support posts. The error In J and e 

P P 

resulting from this omission is difficult to assess accurately. How- 
ever, the physical area of these surfaces was 5% and 2%, respectively, 
of the total interior surface area of the discharge chamber. Thus, one 
might take 7% to be the maximum systematic error in J p resulting from 
the neglected ion currents. Because of these neglected currents, the 



Figure B-1. Ion Source Instrumentation Schematic 
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true value of J p should be larger than the measured value. This implies 
that the true values of Ep and e* should be smaller than the measured 
values according to Eq. 2. A 7% error in J p can result in approximately 
a 13% error in e* . It should be noted, however, that Reference 43 
through 46 indicate that the effective anode area for ion collection at 
a magnetic field cusp should be less than the physical area. Conse- 
quently, it is believed that the measured values of Ep and ej are no 
more than 10% greater than the true values. In addition, the true 
value of f B would be expected to be slightly smaller than the measured 
value as a result of the neglected ion currents. 

Secondary electrons, emitted as a result of ions striking the 

negatively biased surfaces, produce an error in the measurement of the 

ion current to those surfaces resulting in measured values of J that 

c 

are larger than the actual values. The largest secondary electron 
yields at low energies, for the ions and metal surfaces of concern in 
this study, appear to be for argon ions incident on a clean molybdenum 
surface [41]. The secondary electron yield for low energy (< 100 eV) 
argon ions striking molybdenum is approximately 12% [41]. If half of 
the ions produced strike clean molybdenum surfaces then the measured 
value of J p would be approximately 6% greater than the true value as a 
result of secondary electron emission. 

It has been observed, however, that the secondary electron yield 
of a clean surface decreases after exposure to air [59,60]. Reference 
60 concludes that secondary electron yields are less than 1% for low 
intensity and low energy ion incident on a metal electrode which has 
been exposed to air and the ambient gas pressure is in the mtorr range 


A 




- . *1' 




during operation. In the course of the experiments in this investiga- 
tion the discharge chamber surfaces were frequently exposed to air 
between tests. It is believed that the error in the measurement of J p 
resulting from secondary electron emission is less than 1%. 

In the calculations of e p , ej and f g it is assumed that only singly 
charged ions exist. On the basis of this assumption, a doubly charged 
ion leaving the plasma is interpreted as two singly charged ions. The 
effect, however, of these doubly charged ions on e p , ej and f g should 
not be substantial for two reasons. First, the fraction of doubly- to- 
singly charged ions is generally small. Second, the energy required to 
produce one doubly charged ion should not be vastly different than the 
energy required to produce two singly charged ions. The presence of 
doubly charged ions does have a significant effect on the value of the 
neutral density parameter, however, and this problem is addressed later 

in this appendix. 

Of the three systematic errors mentioned above only the neglect of 
the ion currents to anode surfaces and to the cathode support posts 
appears to be significant. The measured values of e p and e* are, there- 
fore, believed to be no more than 10% greater than the true values as a 
result of systematic measurement errors. 

The following analysis was performed to estimate the uncertainty 
associated with the experimentally measured values of e p . The uncer- 
tainty in e p results from the uncertainty in the measurements of the 
independent variables which appear on the right-hand-side of Eq. 2. In 
general, the uncertainty of a quantity y which is a function of the 
measurable independent variables x-j , X 2 » • • * n 9^ ven by [61], 
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where ax., is the uncertainty of the i independent variable. For 
the plasma ion energy cost equation given by*, 

, - (J P"V V D 

P J P 

Eq. B-l becomes, 

Ae p = [( aJ d) + (i^ Aj p) + (3 V^ aV d) . * (B-3) 

Carrying out the partial differentiations, Eq. B-3 becomes, 

4E P = |(^ J P) 2+ (^S) 2+ S^-’) 4 V D ] 2 P . (B-4) 

To use Eq. B-4 the uncertainties in the measurements of the discharge 
current (aJ q ) , the discharge voltage (a V q ) and the total ion current 
produced (AJp) must be determined. These uncertainties result from the 
uncertainties of the digital meters used to make the measurements and 
any variation in the thruster operating set point which may occur while 
the data is being recorded. 

Simpson digital panel meters were used to measure the discharge 
current, discharge voltage and beam current. These meters have an 
uncertainty of 0.1% plus 1 digit. The ion current to the negatively 
biased surfaces was measured using a Kiethly digital multimeter. This 
meter has an uncertainty of 0.5% plus 1 digit. The variation In the 
thruster operating set point was taken to be 1 digit on each of the 


* This form of Eq. 2 is appropriate for the case where J . may 
neglected and the impingement current is included in jJ. 
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digital meters. In addition, 3v was added to the uncertainty In the 
discharge voltage because of the voltage drop across the thermonlc 
cathode wires. With these uncertainties, the uncertainty In e p was 
calculated using Eq. B-4. The results of these calculations are shown 
In Fig. B-2 where the vertical error bars represent the measurement 
uncertainty. The horizontal error bars represent the uncertainty In 
the neutral density parameter and Is discussed In the next section. 

The data In Fig. B-2 Is the same as that In Fig. 6a. The uncertainty 
in the measured values of e p range from approximately 7 to 11%. 


Measurement of m(l-n u ) 


Three systematic errors exist In the measurement of the neutral 
density parameter, m(l-n u ). The first is the effect of doubly charged 
ions on the propellant utilization efficiency. The second Is the use 
of a gas flow meter with argon, krypton and xenon gases that has only 
been calibrated for air. The third is the backflow of neutral atoms 
from the vacuum chamber into the discharge chamber through the grids. 

The presence of doubly charged Ions in the beam leads to artific- 
ally high measurements of the propellant utilization and, consequently, 
low values of the neutral density parameter. Measurements of the 
doubly- to-slngly charged Ion beam current Jg + /Jg were made along the 
thruster centerline. This information can be used to correct the 
propellant utilization according to the equation. 


^ n u^ corrected “ 


1+ \ J B +/J B 


Since the value of Jg /Jg at the centerline is generally higher than 


unmeasured lA .++..+ 
1+ J B ' J B 


(B-5) 
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Figure B-2. Uncertainty in the Plasma Ion Energy Cost Measurements 
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the average value Its use In Eq. B-5 should tend to over correct the 
propellant utilization slightly. The resulting error In the neutral 
density parameter Is believed to be small, however. 

A Teledyne Hastlngs-Raydst flow meter with a range of 0 to 50 
seem was used to measure the propellant flow rate Into the thruster. 

The flow meter was calibrated for air, requiring the output readings to 
be analytically corrected for operation with the gases argon, krypton 
and xenon. The use of the analytical corrections rather than recali- 
brating the flow meter for the different gases Is expected to Introduce 
only a very slight error In the flow rate measurement. Thus, the 
systematic errors In the measurement of the neutral density parameter 
are believed to be negligible. 

The backflow of neutral atoms from the vacuum facility Into the 
discharge chamber was calculated based on measurements of the facility 
pressure made Mm downstream of the thruster. These calculations were 
used to correct the total propellant flow rate Into the discharge 

chamber. 

The accuracy of the flow meter is 1% of full scale. This turns 
out to be the major uncertainty In the determination of the neutral 
density parameter. The uncertainty In the neutral density parameter 
is indicated by the horizontal error bars In Fig. B-2. 

Plasma Property Measurement 

Plasma property measurements were made using one Langmuir probe 
positioned along the discharge chamber centerline and a second probe 
positioned on a magnetic field cusp. Details of the Langmuir probe 
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circuitry used are given in Reference 62. There are two major sources 
of error in using Langmuir probes to obtain plasma property data. The 
first of these deals with obtaining the probe trace itself (i.e., 
measurement of collected current versus probe voltage). The second 
source of error arises from the analysis of the probe trace to obtain 
the desired plasma property information. 

Errors associated with obtaining the probe trace include such 
things as: the variation in work function over the probe surface area, 
secondary electron emission from the probe, probe insulator contamina- 
tion, noise in the plasma, electrical noise in the probe circuitry, 
plasma perturbation by the probe, and magnetic field effects. Plasma 
perturbation by the probe is minimized by using a probe with a small 
surface area. Biasing the probe negative of cathode potential when 
data is not being collected tends to sputter clean the surface and 
minimizes the variation in work function over the probe surface area. 
Changing the probe support insulator frequently every 6 hours) mini- 
mizes the problem of insulator contamination. Noise in the plasma was 
minimized by using a D.C. current heated therminoic cathode. Magnetic 
field effects were minimized by placing the probe on the discharge 
chamber centerline where the magnetic flux density was less than 0.001 
tesla. For the probe positioned on the magnetic field cusp, where the 
magnetic flux density is approximately 0.2 tesla, the magnetic field 
effects are substantial. In this case, only the electron temperature 
and primary electron energy data obtained from these traces are believed 
to be meaningful . 


The probe traces were analyzed using a non-linear numerical curve 
fitting routine similar to the one described by Beattie [6]. This 
routine assumes an electron population that is characterized by a Max- 
wellian plus mono-energetic (primary) electron energy distribution. 

Probe traces were digitized using an HP-7470A graphics plotter together 
with an HP-85 mini -computer. The plotter may be used as a digitizer 
when equipped with a special digitizing sight. Accurate digitizing is 
crucial to obtaining good primary electron information. 

The errors introduced during probe trace digitizing and analysis 
may be determined by generating idealized probe traces with known 
plasma properties and then analyzing these traces in the usual manner. 
Eight such idealized probe traces were generated and analyzed, covering 
a wide range of plasma conditions. The results indicated that the data 
analysis can accurately distinguish primary-to-Maxwellian electron 
density ratios as small as 0.2% provided the Maxwellian electron temper- 
ture is low. The most difficult traces to analyze are those correspond- 
ing to low primary-to-Maxwellian electron density ratios in a plasma 
with a high Maxwellian electron temperature. This particular plasma 
condition, however, was not observed experimentally in the course of 
this investigation. The results in Chapters III, IV and Appendix A of 
this report indicate that when the primary-to-Maxwellian electron 
density ratio Is small the Maxwellian electron temperature is low and 
when the electron temperature is high the primary-to-Maxwellian elec- 
tron density ratio is large. These types of probe traces are accurately 
analyzed with the data reduction system used in this investigation. 

The idealized probe trace analysis indicated that under most conditions 


the following errors may be expected from the digitizing and curve 


fitting procedures: 

Plasma Potentional i 0.3 v 

Maxwellian Electron Density, n M t 5% 

Maxwellian Electron Temperature, t 3% 

Primary Electron Density, n p t 15% 

Primary Electron Energy, Ep i 5% 

V"m : 3* 

The overall uncertainty associated with the plasma property measure- 
ments are believed to be the following: 

Plasma Potential 1 1.0 v 

Maxwellian Electron Density, n^ t 20T 

Maxwellian Electron Temperature, T M t 10% 

Primary Electron Density, n t 30% 

r 

Primary Electron Energy, E t 10% 

r 

y n M t 30% 


Measurement of Doubly Charged Ions 
Measurements of doubly charged ions In the beam were made using an 
E x B probe similar to that described in Reference 42. Errors associ- 
ated with these measurements include: the accuracy of the picoammeter 
used to read the probe current, selection of the proper probe deflection 
voltage to obtain the maximum signal corresponding to the particular 
charged specie of Interest, and variation In the thruster operating 
conditions during data collection. 
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Probe current readings were made using a Keithly model 410A 
plcoanmeter with an accuracy of ± 4% of full scale. Selection of the 
proper probe deflection voltage Is a fairly difficult thing to do 
accurately because of the slow response time of the plcoammeter 
(between 0.4 and 12s). This Is especially true for currents In the 
lowest ranges of the ammeter. Errors resulting from the uncertainty 
in probe deflection voltage and variation in thruster operating condi- 
tions are believed to Introduce a maximum error of 10* In the measured 
values of the single and double Ion currents. The total uncertantles 
in the measured values of the doubly-to-slngly charged Ion current 
ratio, resulting from the sources of error mentioned above, are indi- 
cated in Fig. B-3 by the vertical error bars. The horizontal error 
bars indicate the uncertainty In the corrected propellant utilization 
efficiency as mentioned earlier. 


OOUBLE-TO-SINGLE ION CURRENT RATIO 
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APPENDIX C 
NOMENCLATURE 

A g a Area of grids through which the ion beam Is extracted (m 2 ) 

B = mC 

o o 

C Q ■ Primary electron utilization factor defined by Eq. 26 (A eq . ) ~ 1 

E p ■ Primary Electron Energy (eV) 

e * Electronic Charge (1.6 x 10" 19 coul.) 

f A * Fraction of ion current to anode surfaces 

f B * Extracted ion fraction 

f c = Fraction of 1;> . current to cathode surfaces 

J A = Ion current to anode potential surfaces (A) 

Jg * Ion beam current (A) 

Jg * Singly charged ion beam current (A) 

Jg + * Doubly charged ion beam current (A) 

Jq = Ion current to cathode potential surfaces (A) 

* Discharge current (A) 

J E * Cathode Emission Current (A) 

J* « Total production rate of excited neutral atoms by primary 
ex electrons - expressed as a current (A) 

Jimp * Accelerator grid ion impingement current (A) 

J^ = Production rate of jth excited state expressed as a current (A) 

J L = Primary electron current to the anode (A) 

J M s Maxwellian electron current to the anode (A) 

Jp = Total ion current production rate expressed as a current (A) 


124 


U P 

*T 

p,m 


to, 


M c 

M.. 


M| 

m ( 

m i 

• 

m 

rn 


'Po 


p o + 


p 0+t , 


* current d (A) l0n r “ e by pr1nBry e,ectr0 " s expressed as a 
rate of sin 9'y charged Ions expressed as a 

current (A) 

Production rate of doubly charged ions expressed as a 
current (A) 

= current^(A)^° n rate by Maxwe11ian electrons expressed as a 

s Primary electron containment length (m) 

5 Final spacecraft mass (kg) 

8 Generator mass (kg) 

8 Initial spacecraft mass (kg) 

! Payload mass (kg) 

Initial propellant mass (kg) 

Electron Mass (kg) 

Ion mass (kg 

Thruster propellant flow rate (A eq.) 

Total spacecraft propellant flow rate (kg/s) 

Total electron density (m -3 ) 

Total ion density (m~ 3 ) 

Maxwellian electron density (m“ 3 j 
Primary electron density (m“ 3 ) 

Neutral atom density (nr 3 ) 

Neutral atom loss rate (A eq.) 

Generator Power (W) 

rate f4Ct0r f0r 1on,2at1on ° f "eotral 
neutral Itoms^raVs)* faCt ° r d ° Ub ' e 1o " 1m,on » f 
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P ++ = Primary electron rate factor for the production of double 

+ Ions from single Ions (m 3 /s) 

Q + = Maxwellian electron rate factor for Ionization of neutral 
0 atoms (m 3 /s) 

Qq + = Maxwellian electron rate factor for double Ionization of 
neutral atoms (m 3 /s) 

q+ + = Maxwellian electron rate factor for the production of 
+ double ions from single Ions (m 3 /s) 

T. = Maxwellian electron temperature at anode surface (eV) 

H 

T = Maxwellian electron temperature in bulk plasma (eV) 

M 

t = Mission duration (s) 
u * Ion exhaust velocity (m/s) 

U * Lumped excitation energy (eV) 

CA 

U + ■ Ionization energy (eV) 

Uj * Excitation energy of jth excited state (eV) 

U * Lowest excitation energy (eV) 

= Anode sheath voltage (v) 

V = Plasma potential from which electrons emitted by the cathode 
c are accelerated to become primary electrons 

V D = Discharge (Anode) voltage (v) 

V.. * Net accelerating voltage (v) 

N 

V + = Screen grid supply voltage (v) 

V = Accelerator grid supply voltage (v) 

v, * Bohm velocity (m/s) 

0 

v = Electron velocity (m/s) 
e 

v a Primary electron velocity (m/s) 

P 

v„ = Neutral atom velocity (m/s) 
o 

V- = Volume of ion production region (m 3 ) 
x = Independent variable 
y = Dependent variable 
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Power plant specific mass (kg/w) 

Discharge current uncertainty (A) 

Ion production current uncertainty (A) 

Characteristic velocity (m/s) 

Discharge voltage uncertainty (V) 

Uncertainty of Independent variable 

Uncertainty of dependent variable 

Average beam ion energy cost (eV/beam ion) 

Average energy of Maxwellian electrons leaving the plasma 
at the anode (eV) 

Average plasma Ion energy cost (eV/plasma ion) 

Baseline plasma ion energy cost (eV/plasma Ion) 

Average plasma Ion energy cost considering ionization and 
excitation processes only (eV) 

Thruster electrical efficiency 

Propellant utilization efficiency 

Total excitation collision cross section (m 2 ) 

Total excitation collision cross section at the primary 
electron energy (m 2 ) 

Excitation collision cross section of the jth state (m 2 ) 

Excitation collision cross section of the jth state at the 
primary electron energy (m 2 ) 

Total inelastic collision cross section at the primary 
electron energy (m 2 ) 

Ionization collision cross section (m 2 ) 

Ionization collision cross section at the primary electron 
energy (m 2 ) 

Transparency of the accelerator grid to neutral atoms 
Transparency of the screen grid to ions 




